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SECTION  1 

; INTRODUCTION 

This  was  a program  directed  toward  an  improved  understanding  of  the 
sources  and  energization,  transport  and  loss  processes  associated 
with  the  low-energy  (keV)  plasmas  in  the  earth's  magnetosphere, 
particularly  as  they  relate  to  ionospheric  perturbations  and 
ionospheric-magnetospheric  coupling.  ^ Under  this  program^ the  Lockheed 
Palo  Alto  Research  Laboratory  (LPARL)*>  conducted  a low-energy  particle 
experiment  on  the  S3-3  spacecraft,  ^which  was  launched  in  the 

summer  of  1976s-*  The  experiment  is  still  operating  successfully  and 
has  made  several  major  discoveries  about  previously  unknown  ionospheric 
acceleration  mechanisms  and  the  composition  of  the  ring  current. 
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Section  2 

BACKGROUND/ RELEVANCE 

A number  of  Navy  systems  either  depend  upon,  or  are  affected  by,  the  iono- 
sphere and  the  magnetosphere.  Communication  links  utilize  ionospheric 
propagation  paths  which  are  significantly  disturbed  by  the  precipitation 
of  energetic  magnetospheric  plasma  in  the  keV  range,  particularly  in  the 
polar  regions.  The  heating  and  expansion  of  the  upper  atmosphere  due  to 
this  precipitation  can  increase  drag  and  affect  the  orbits  of  low  altitude 
navigational  satellites.  The  interaction  of  the  hot  magnetospheric 
plasma  with  synchronous  satellites  leads  to  a variety  of  pheonomena  associated 
with  spacecraft  charging  which  degrades  the  performance  of  communications 
satellites  in  operational  systems. 

A new  dimension  has  recently  been  added  to  our  ability  to  study  the  above 
phenomena  with  the  development  of  energetic  ion  mass  spectrometers  which 
for  the  first  time  allow  for  the  identification  of  the  ion  species  in  the 
measured • plasma . This  capability  has  produced  some  remarkable  results, 
principal  of  which  was  the  discovery  that  0+  is  on  occasion  the  dominant 
species  in  the  few  keV  range.  Since  the  0+  ions  have  a much  smaller  range 
than  protons  of  the  same  energy,  they  have  dramatically  different  interac- 
tion characteristics  with  both  the  upper  atmosphere  and  with  sensitive 
spacecraft  surfaces  and  therefore  the  capability  to  identify  the  ion  spe- 
cies is  required  for  an  understanding  of  the  plasma  effects.  The  ONR-118 
spectrometer  provides  this  identification  via  the  utilization  of  both 
velocity  and  energy  analysis  on  the  measured  ions.  It  covers  the  energy 
range  0.6  to  16  keV  per  unit  charge  for  ions  and  0.1  to  20  keV  for  electrons. 


2-1 


LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 

itCINlll  » » X % I I • % ft  ! M ( I C O M M N T INC 
A Of  lOriNIIO  I l I f I I M (OIIOimoN 


LSMC/D673078 


Section  3 
RESULTS 

Several  important  discoveries  have  resulted  from  the  experiment.  The 
principal  result  is  the  discovery  of  large  fluxes  of  energized  ionospheric 
ions  streaming  up  the  magnetic  field  lines  from  the  ionosphere.  This 
discovery  establishes  the  existence  of  a major  new  ionospheric-magnetos- 
pheric  coupling  phenomenon.  The  characteristics  of  the  phenomenon  point 
to  parallel  electric  fields  as  the  immediate  cause  of  the  accelerated 
ions.  These  results  and  those  of  the  ONR  companion  experiment  on  S3-3, 
which  makes  in  situ  measurements  of  the  electric  fields,  put  on  a firm  basis 
the  many  recent  speculations  and  inferences  relating  to  large  electric 
potential  differences  parallel  to  the  geomagnetic  field  lines  in  the  magnet- 
osphere. 

Another  important  result  was  the  discovery  of  a second  type  of  ionospheric 
acceleration  mechanism  which  preferentially  accelerates  the  perpendicular 
component  of  the  ion  velocity. 

A third  important  result  was  the  first  direct  measurement  of  the  composi- 
tion of  the  storm-time  ring  current.  In  the  energy  range  from  1-16  keV 
the  main  phase  ring  current  has  been  found  to  be  primarily  composed  of  H+ 
and  0+  ions . 

f 

Some  of  our  more  recent  results  derive  from  detailed  case  studies  and  stat- 
istical studies  of  the  phenomena  which  were  discovered  in  the  initial  phases 
of  the  experiment.  A quantitative  analyses  of  the  signatures  in  the  ion 
and  electron  distribution  functions  is  leading  to  detailed  information  on 
the  geometry  and  physical  processes  involved  in  the  parallel  electric  field 
regions  (see  Appendix  A and  B) . A statistical  study  of  downward  flowing 
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ion  events  has  led  to  the  conclusion  that  the  ionosphere  Is  a principle 
contributor  to  the  circulating  ion  clouds  which  dominate  the  plasma  en- 
vironment in  the  vicinity  of  the  geostationery  orbit. 


These  results  have  been  described  in  a series  of  publications  and  pre- 
sentations which  are  included  in  the  following  bibliography.  Preprints 
and  reprints  of  the  most  recent  articles  are  included  in  the  Appendix. 
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ABSTRACT 


Measurements  of  energetic  electrons  (0.07  < E < 24  keV) 

and  ions  (0.3  <_  E .1  16  keV)  on  the  S3-3  satellite  show  features 

which  are  interpreted  in  terms  of  parallel  potential  drops  of 

up  to  several  kV  magnitude  existing  simultaneously  both  above 

and  below  the  satellite.  This  leads  to  the  inference  that  the 

satellite  is  within  an  auroral  acceleration  region.  Two  events 

of  this  type  are  examined  in  which  the  satellite  is  in  the  7000- 

8000  km  altitude  range  and  at  auroral  latitudes.  Both  events  are 

interpreted  as  traversals  of  an  "inverted  V"  acceleration  region 

of  broad  latitudinal  width  (2°  - 4°).  In  both  cases,  one  can 

3 

infer  a vertical  extent  to  the  acceleration  region  of  > 10  kn. 
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I.  INTRODUCTION 


V 


Previous  results  have  suggested  that  at  least  some  types  of  auroral  electron 

acceleration  processes  occur  within  large  scale  Ionospheric  structures  contain- 

4 

ing  potential  drops  of  the  order  of  10  volts  aligned  parallel  to  the  magnetic 
field  and  situated  at  altitudes  above  those  which  have  been  repeatedly  sampled 
by  experiments  on  rockets  and  low  altitude,  polar-orbiting  satellites  [Gurnett, 
1972;  Evans,  1974;  Swift,  1975;  Kauffman  et  al.,  1976,  Falthammer,  1977]. 

The  S3-3  satellite  (8,000  km  apogee,  240  km  perigee,  97.5"  inclination) 
is  now  sampling  this  region  of  the  magnetosphere  and  has  discovered  that  such 
structures  arc  indeed  a commonly  occurring  phenomena  at  altitudes  of  about 
1 Rj.  at  auroral  latitudes.  Some  of  the  manifestations  of  these  phenomena 
are: 

1.  Upward  flowing,  narrowly  collimated  beams  of  energetic  (keV) 
field  aligned  ions,  indicative  of  electrostatic  acceleration 

of  ambient  ionospheric  ions.  [Shelley  et  al.,  1976,  Ghiclmctti 
et  al.,  1978] 

2.  Occasional  conical  pitch  angle  distributions  in  the  energetic 
ions  indicating  additional  processes  providing  a perpendicular 
as  well  as  a parallel  component  to  the  ion  accelerat ion.  [Sharp 
et  al.,  1977;  Ghiclmctti  et  al.,  1978]. 

3.  Signatures  of  keV  potential  drops  both  above  and  below  the 
satellite  in  the  energetic  electron  populations  [Mizcra  and 
Fennell,  1977;  Cladis  and  Sharp,  1977a,  b]. 

4.  Narrow  spatial  regions  with  measured  d.c.  electric  fields  > 

2 

10  mV/m  which  are  referred  to  as  electrostatic  shocks  by  the 
ers 

Berkeley  experiment/  [Mozer  et  al.,  1977;  Torbert  and  Hozer,  1978]. 
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IV.  SUMMARY  AND  CONCLUSIONS 
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In  this  paper  we  show  some  examples  of  the  plasma  observations  from 
within  these  ionospheric  structures,  focusing  primarily  on  the  use  of  the 
electron  pitch  angle  distributions  as  a diagnostic  tool  to  probe  for  indica- 
tors of  their  geometry,  particularly  their  vertical  extent.  A detailed  dis- 
cussion of  the  use  of  this  type  of  analysis  on  another  event  has  been  pre- 
sented by  Cladis  et  al.  (1977). 

The  Lockheed  experiment  on  S3-3  consists  of  three  ion-mass-spcct rometers 
and  four  magnetic  electron  spectrometers  which  have  been  described  in  previous 
publications  (Shelley  et  al.,1976,  sharp  ot  al.,  1977  ].  The  satellite 

is  spinning  at  approximately  3 rpm  with  its  spin  axis  perpendicular  to  the 
orbital  plane  and  the  instruments  are  mounted  with  their  view  directions 
perpendicular  to  the  spin  axis,  providing  almost  complete  pitch  angle  scans 
of  the  measured  fluxes.  The  angular  acceptance  range  for  the  spectrometers 
is  approximately  6®  full  width.  The  electron  spectrometers  have  broad  energy 
bands  and  nearly  rectangular  response  functions  so  as  to  include  almost  all 
electrons  with  .07  <_  E <_  24  kcV  in  the  four  measured  intervals  (See  Table  1). 
They  are  sampled  twice  per  second  providing  approximately  9®  angular 
resolution.  The  ion  spectrometers  each  have  4 energy-per-charge  settings 
which  arc  stepped  every  16  seconds.  These  are  also  listed  in  Table  1.  The 
ion  spectrometers  provide  a complete  truss -pet -charge  sweep  at  each  of  three 
energy-per-charge  values  every  1 second. 

Since  the  satellite  spins  at  a.  3 rpm  and  it  requires  several  spins  to 

be  able  to  differentiate  characteristic  pitch  angle  structure  from  spatial 

or  temporal  variations,  we  focus  these  studies  on  ionospheric  structures  with 

■ broad  latitudinal  extent  which,  as  wc  will  sec,  may  have  relatively 

weak  parallel  electric  fields  (of  the  order  of  millivolts  per  meter)  ertending 

3 

over  vertical  dimensions  of  >\.  10  km.  Thoso  arc  possibly  a different  class  of 

phenomena  than  those  being  studied  by  the  Berkeley  group  [Mozoc  et  al.,  1977] 

whose  experiment  more  naturally  focuses  on  high  electric  field  phenomena 
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f'lO*’  mv/m)  which  arc  found  to  occur  over  narrow  few  km)  latitudin.il  di.v.en- 
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•lons[Torbert  and  Mozer  1978].  The  two  claascs  of  phenomena  are  generally 
associated  spatially  but  arc  not  identical.  An 

example  of  this  spatial  association  is  given  in  Figure  1 which  shows  data 
obtained  over  the  northern  hemisphere  auroral  zone  on  29  July  1970.  The 
abscissa  shows  universal  time,  altitude  in  kilometers,  invariant  latitude, 

Magnetic  local  time,  and  L value.  The  four  lowest  panels  show  the  logarithm 
of  the  counts  per  half-second  counting  interval  for  the  electron  spectrom- 
eters. The  panel  labeled  PITCH  shows  the  pitch  angle  of  the  look  direction 
of  the  spectrometer.  The  next  four  panels  show  the  logarithm  of  the  counts 
from  ions  with  m/q  = 1,2,4,  and  16,  respectively,  summed  once  per  second 
from  selected  output  channels  from  all  three  of  the  mass  spectrometers  giv- 
ing an  approximate  measure  of  the  relative  flux  of  the  relevant  species.  The 
top  panel  shows  the  energy  step  of  each  of  the  three  mass  spectrometers.  The 
spin  numbers  arc  labeled  for  reference.  For  example,  for  the  16-second  period 
associated  with  spin  12  the  three  spectrometers  were  set  at  E/q  values  of 

I. 28,  4.5,  and  16.0  keV,  respectively,  and  the  m/q  = 1 plot  shows  the  sum  of 
the  proton  counts  at  these  three  energies.  In  fact,  on  this  spin  the  observed 
counts  derived  almost  entirely  from  the  high  energy  spectrometer  indicating  a 
peak  in  the  spectrum  at  E >,  16  keV. 

In  examining  Figure  1 one  sees  two  regions,  labeled  Region  I and  Region 

II,  each  of  which  contains  particle  fluxes  with  common  signatures  extending 
over  a wide  latitudinal  range.  Region  I shows  symmetric,  90®  peaked  pitch 
angle  distributions  in  the  three  lowest  energy  electron  channels  which,  as 
will  be  shown,  are  indicative  of  a broad  region  of  parallel  electric  fields 
located  above  the  spacecraft.  Region  II  is  a broad  region  of  upward  flowing 
Ions  associated  with  a drop-out  in  the  low  energy  electron  fluxes  which  is 
indicative  of  a potential  drop  below  the  spacecraft. 

Both  of  these  regions  are  contained  within  a larger  scale  region  in 
which  the  electron  fluxes,  as  measured  by  the  Aerospace  Corporation's  electro- 
static analyzers,  show  the  characteristic  signature  of  an  "inverted  V"  event 

A-5 


LMSC/D673078 


References 


LMSC/D673078 


when  plotted  on  a grey  scale  energy-time  spectrogram.  [Mizera  and  Fennell, 
1978;  the  spectrogram  for  this  event  is  shown  in  Figure  3 of  this  reference.] 
Within  this  large  scale  "inverted  V"  region  four  discrete,  narrow,  electro- 
static shocks  have  been  reported  by  the  Berkeley  group  [Mozcr  ct  al.,  1977] 
at  the  locations  indicated  on  the  lowest  horizontal  axis  in  Figure  1.  Upward 
flowing  energetic  ions  are  occasionally  observed  at  the  locations  of  the  elec- 
trostatic shocks  [Mizera  ct  al.,  1977]  when  the  plasma  analyzers  are  fortui- 
tously oriented  downward  at  the  times  of  passage  through  the  narrow  structures 
On  this  occasion  the  Lockheed  spectrometers  saw  no  upward  flowing  ions  until 
the  satellite  reached  Region  II. 

As  indicated  above,  this  report  will  focus  on  the  broad  scale  structures. 
Preliminary  indications  are  that  these  regions  are  generally  spatially  asso- 
ciated with  enhanced  low  frequency  noise  or  turbulence  in  the  electric  field 
data  [Torbert  and  Mozer,  1978].  The  signatures  of  the  parallel  potential 
drops  within  the  structures  consist  of  the  upward  flowing  ion  beams  plus  cer- 
tain features  of  the  electron  pitch  angle  distributions  to  be  described  below. 
Signatures  of  potential  drops  above  and  below  the  satellite  are  often  observed 
simultaneously  indicating  a substantial  vertical  extent  to  the  structures. 

As  discussed  by  Chiclmctti  et  al.  [1978]  the  ion  beams  often  show  the  com- 
bined effects  of  perpendicular  and  parallel  accelerations  as  well  as  the 
transfer  of  energy  between  the  two  components  by  the  mirroring  action  of  the 
magnetic  field.  This  complication  plus  the  time  lags  in  the  ion  data  due  to 
the  substantial  times  of  flight  involved  (e.g.,  1 kcV  0+  ions  require  over 
9 secs  to  travel  1000  km)  make  the  ion  data  less  useful  for  quantitative 
analysis  in  possibly  time-varying  events  than  the  rapidly  moving  electrons. 
Therefore,  in  these  initial  studies  we  will  focus  on  the  electron  signatures. 
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For  the  following  discussion  we  define  primary  electrons  as  those  electrons 

entering  the  acceleration  region  from  above  and  those  electrons  which  have 

passed  through  the  acceleration  region  and  are  reflected  below  the  spacecraft 

before  interacting  with  the  atmosphere.  We  define  degraded  primary  electrons 

as  those  primary  electrons  which  have  lost  energy  in  the  atmosphere  below 

the  acceleration  region,  and  we  define  secondary  electrons  as  those  electrons 

which  have  been  produced  in  the  atmosphere  below  the  acceleration  region  by 

the  primary  electrons.  As  illustrated  in  Figure  2,  the  electrons  observed  in 

the  presence  of  a parallel  upward-directed  electric  field  above  the  spacecraft 

fall  into  two  general  categories;  primaries  and  degraded  primaries  with  energy 

greater  than  the  potential  drop  above  the  spacecraft,  (II  > , and  secondaries 

and  degraded  primaries  with  energy  less  than  the  potential  drop  (E  < $ ).  The 

class  is  characterized  by  an  up-down  symmetry  in  their  pitch  angle  distributions 

since  the  upcoming  electrons  from  the  atmosphere  are  all  reflected  by  * . 

S«ince  the  primary  electrons  are  all  accelerated  to  energies  greater  than  ,*  , a 

A 

detector  whose  energy  band  is  entirely  below  $ will  shown  an  up-down  svr.-.otry 
as  indicated  in  Figure  2a.  In  addition  to  the  short-lived  particles  topologi- 
cally connected  to  the  atmosphere  (shaded  region  in  Figure  2a),  one  often  sees 
a quasi-stable  component  of  the  fluxes  in  this  energy  range  which  are  trapped 
between  the  magnetic  mirror  below  the  spacecraft  and  the  electric  mirror  above 
(angular  regions  indicated  by  circles  in  Figure  2a).  This  population  is  fed 
by  scattering  and  by  fluctuations  in  the  electric  fields  and,  because  of  inten- 
sity fluctuations  in  the  flux  levels  of  the  primaries,  can  build  up  to  levels 
greater  than  the  Instantaneous  levels  observed  in  the  loss  cone.  Examples  of 
this  class  of  signatures  are  the  Region  I electron  fluxes  in  CMEA,  CMKB,  and 
CMEC  shown  in  Figure  1 and  the  CMEA  fluxes  in  the  central  portion  of  the  event 
illustrated  in  Figure  7. 
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The  electrons  with  E > *A»  generally  exhibit  a well-defined  loss  cone 
which  destroys  the  up-down  symmetry  and  establishes  their  identification. 

The  loss  cone  is  both  widened  and  deepened  by  the  potential  below  the  space- 
craft [Cladis  et  al.,  1977]  with  the  widening  simply  related  to  the  magnitude 
of  this  potential  through  the  expression: 


. 2 

sin  a. 


ES  + *B 


(1) 


where  a is  the  angular  half  width  of  the  loss  cone.  6 is  the  kinetic 

energy  gained  in  the  parallel  field  below  the  satellite,  B and  B are  the 

S T 

magnitudes  of  the  magnetic  field  at  the  location  of  the  satellite  and  the 
top  of  the  atmosphere  respectively,  and  E is  the  measured  energy  of  the 
electrons  at  the  location  of  the  satellite. 

The  primary  electrons  with  energies  only  moderately  greater  than  i 

A 

also  exhibit  an  angular  cutoff  in  the  vicinity  of  90°.  This  can  be  inter- 
preted as  a result  of  the  action  of  the  first  adiabatic  invariant: 


B„ 


Es  - 


“"2«  C • ^ ES 


(2) 


where  c»c  is  the  angular  location  of  the  cutoff,  A is  the  kinetic  energy 

gained  in  the  parallel  field  between  the  satellite  altitude  and  (see 

Figure  2)  and  B^  is  the  magnetic  field  strength  at  Z . This  equation  has  two 

quantities  which  are  generally  unknowns,  A and  B. . B,  is  of  particular 

interest  since  it  is  related  to  the  vertical  scale  size  of  the  ionospheric 

structure  (labeled  d in  Figure  2)  which  can  serve  to  help  differentiate 

between  various  models  for  the  cause  of  the  parallel  potential  drops  (e.g. , 

* 

double  layers,  oblique  shocks,  anomalous  resistivity)  which  have  character- 
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istically  different  scale  sizes.  In  this  work  we  will  assume  isotropy  of 
the  Incident  spectrum  (for  Z > Z^)  and  utilize  simultaneous  measurements  of 
electrons  in  the  different  energy  ranges  to  provide  information  on  both  un- 
knowns in  Equation  (2)  thereby  inferring  values  for  the  parameter  d , as 
well  as  obtaining  estimates  of  ^ and  from  the  measured  cutoffs  in  the 
pitch  angle  distributions  at  and  a^. 

Since  we  are  utilizing  broad  band  electron  detectors  some  uncertainty 

arises  in  the  assignment  of  E in  these  expressions  for  a and  Q . To 

*>  CL 

improve  the  accuracy  of  the  estimates. of  $ , a , and  d one  can  fit  an 

A D 

assumed  distribution  function  to  the  detector  response  or  otherwise  use 

supplementary  information  from  other  detectors.  We  have  taken  this  approach 

to  estimate  the  energy  E to  be  associated  with  the  pitch  angle  q correspond- 

5 C 

ing  to  50Z  of  the  total  reduction  in  the  count  rate  in  the  region  of  the 
angular  cutoff. 

A more  general  and  extended  discussion  of  these  and  other  signatures 
of  parallel  electric  fields  in  the  particle  distribution  functions  lias  been 
presented  by  Kauffman  ct  al.  [1976],  and  Whipple  [1977]. 

II.  JULY  29,  1976,  11:28  UT  EVENT 

As  discussed  above,  the  up-down  symmetry  in  the  pitch  angle  distribu- 
tions in  CMEA,  CMEB,  and  CMEC  in  Region  I indicates  that  they  arc  respond- 
ing primarily  to  electrons  with  energies  < a . This  implies  that  the  pri- 
■ary  electrons  are  accelerated  out  of  the  detector's  energy  range  and  estab- 
lishes a lower  limit  of  $ i 5 keV  in  this  region.  CMED  which  is  respond- 

A 


ing  to  the  primary  portion  of  the  electron  spectrum  shows  the  90*  mininuims 

in  the  pitch  angle  distributions  expected  to  result  from  $ as  discussed 

A \ 

in  Section  I.  The  transition  from  90*  maximum  to  90*  minimum  type  pitch 
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angle  distributions  evidently  occurs  in  the  lower  end  of  the  CMKD  energy 
window  since  the  observed  90*  minlmums  are  occasionally  somewhat  obscured 
by  electrons  which  retain  the  90*  peaked  angular  distributions  characteristic 
of  the  lower  energy  channels  (see  Figure  1).  The  CMKD  pitch  angle  distribu- 
tion for  spin  3+  (i.o.,  the  second  of  the  two  pitch  angle  scans  on  spin  3) 
is  illustrated  in  expanded  form  as  the  upper  curve  of  Figure  3.  The  average 
value  of  a as  defined  above  for  this  distribution  is  63*.  In  order  to  at- 

v 

tempt  to  determine  the  uncertainty  introduced  by  the  contribution  from  the 
90*  peaked  low  energy  electrons,  we  note  from  Figure  1 that  the  shape  of 
the  90*  maxima  in  the  angular  distributions  is  generally  not  a strong  func- 
tion of  energy.  This  is  illustrated  in  more  detail  for  spin  3+  in  Figure  4. 
We  therefore  subtract  a quantity  with  the  shape  of  the  measured  CMKC  pitch 
angle  distribution  from  the  CMKD  distribution  in  Figure  3 in  order  to  evalu- 
ate the  uncertainty  such  a contribution  might  introduce  to  the  estimated 
value  of  a . The  maximum  magnitude  such  a contribution  could  have  is  that 

Vo 

required  to  reduce  the  measured  CMKD  90°  flux  to  zero.  The  residual  CMKD 
flux  after  the  subtraction  of  such  a distribution  is  illustrated  by  the 
lowest  curve  in  Figure  3.  A similar  curve  corrected  for  a 90*  peaked 
contribution  of  half  this  magnitude  is  illustrated  by  the  Center  curve. 

Both  of  these  distributions  yield  average  values  of  63'1  equivalent 
to  that  derived  from  the  uncorrected  (solid)  curve.  Thus  the  measured  value 
of  this  quantity  is  insensitive  in  this  case  to  the  magnitude  of  the  correc- 
tion and  can  be  used  to  estimate  from  Equation  (2)  and  the  parameter  d 
••  discussed  in  Section  I. 

Within  the  range  of  the  CMF.D  detector  the  electron  flux  during  this 
period  is  generally  decreasing  with  increasing  energy  [P.  F.  Mizcra,  private 
coomunicatlon]  and  the  midpoint  of  the  energy  band  (15.4  keV)  can  be  taken 
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as  an  upper  limit  on  E . Using  the  lower  limit  of  5 keV  on  a derived  from 

s TA 

the  CMEC  pitch  angle  distribution,  we  obtain  Bc/B^ ^ 1.2  and  d j 800  km 
assuming  a 1/r^  magnetic  field  dependence. 

An  alternative  lower  limit  can  be  obtained  provided  there  is  some 
residual  primary  flux  at  90°  within  the  bandpass  of  the  detector  which  has  not 
been  excluded  by  the  action  of  $ (for  example  one  of  the  two  upper  curves 
in  Figure  3).  In  these  circumstances,  electrons  with  energies  at  least  as  low 
as  the  upper  end  of  the  CMED  energy  window  (23.5  keV)  can  reach q > 90°  at 
the  satellite  altitude.  Thus  in  Equation  (2)  we  can  set  = 90®  and  Eg  < 23.5 
keV  for  a limit  on  Bg/B^  i 1.3  and  d j 1100  km. 

The  parallel  potential  drop  below  the  spacecraft  in  region  I is  of 
relatively  low  magnitude.  The  width  of  the  CMED  loss  cone  (c^  = 19°)  shown 
for  spin  3+  in  Figure  3 is  equal  to  its  expected  value  in  the  absence  of  such 


a potential.  No  upward-flowing  ions  (with  E > 0.5  keV)  were  observed  in  the 
Lockheed  mass  spectrometers.  The  Aerospace  electrostatic  analysers  did  in 
fact  see  weak  upstreaming  ion  fluxes  in  the  100  cV  range.  [Mizera  and 
Fennell,  1978]  so  a weak  potential  probably  did  exist  below  the  spacecraft, 
but  its  magnitude  was  too  low  to  be  detected  with  the  techniques  utilized 
here. 

Region  II  on  the  other  hand  is  characterized  by  the  signatures  of  a 
large  As  discussed  in  Section  I,  one  sees  continuous  upstreaming  ion 

beams  with  energies  up  to  16  keV.  The  electron  loss  cones  are  also  signifi- 
cantly widened  with  respect  to  Region  I (see  Figure  1).  This  is  illustrated 
for  spin  7+  in  more  detail  in  Figure  5 where  for  CMED  ■ 28®  corresponding 
to  a ^ : 13  keV. 

Another  indicator  of  the  strong  in  Region  II  is  the  dramatic  suppres- 
sion of  the  low  energy  secondary  and  degraded  primary  electrons  in  this 
region  relative  to  Region  I.  A potential  drop  above  the  spacecraft  is  also 
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implied  by  this  signature  since  the  primary  spectrum  is  accelerated  to  higher 
energies,  but  in  general  this  is  not  sufficient  to  suppress  the  low  energy 
electron  fluxes  since,  depending  on  the  primary  spectrum,  the  fluxes  of  up- 
coming and  reflected  secondaries  in  a given  low  energy  region  can  exceed  the 

fluxes  of  primaries  that  existed  in  that  region  prior  to  the  onset  of  the  * . 

A 

The  *fl  on  the  other  hand  acts  both  to  reflect  the  lowest  energy  secondaries 

and  reduce  the  level  of  the  degraded  primaries  by  the  factor  E /(E  + $ ) 

S S B 

relating  the  flux  at  the  top  of  the  atmosphere  to  the  flux  at  the  satellite 
which  enters  through  Liouville's  theorem. 

In  contrast  to  the  potential  below  the  spacecraft,  the  potential  above 
the  spacecraft  has  not  changed  substantially  between  Region  I and  Region  II. 
Even  though  the  magnitude  of  the  low  energy  secondaries  and  degraded  primary 
electrons  is  suppressed,  when  detectable  (i.e.,  in  CMEC)  they  still  exhibit 
the  up-down  symmetry  characteristic  of  $ :>  5 keV  (see  Figure  6).  (Later, 

toward  the  end  of  Region  II  [“'•spin  12]  one  secs  in  CMEC  the  transition  to  a 
loss  conc/ant iloss  cone  asymmetry  indicating  that  $ is  decreasing  to  the 
OlEC  energy  range  [see  Figure  1].)  The  CMED  distribution  for  spin  7+ 
illustrated  in  Figure  5 shows  a 90®  minimum  with  a = 62°,  about  the  same 
as  in  spin  3+,  again  approximately  independent  of  a correction  for  a low 
energy  component  peaked  at  90®.  From  the  Aerospace  spectrogram  (op.  cit.) 
one  finds  that  the  energy  spectrum  of  the  primary  electrons  in  this  period 
has  not  deviated  appreciably  from  Region  I.  Thus  we  can  conclude  that  the 
electric  field  geometry  above  the  spacecraft  is  approximately  the  same  as  in 
Region  I ( d extending  to  £ 10^  km)  but  that  the  field  below  the  spacecraft 
has  dramatically  increased. 
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III.  SEPTEMBER  15,  1976,  1058  UT  EVENT 

The  July  29  event  described  in  Section  II  was  selected  primarily  to  illus 
trate  the  relationship  of  the  broad  parallel  electric  field  regions  to  the 
electrostatic  shock,  regions  observed  by  Mozer  et  al . (1977).  It  was  a rela- 
tively complicated  event  with  overall  potential  drops  of  larger  magnitude 
than  are  typical.  In  this  section  we  will  describe  a simpler,  less  intense, 
and  isolated  event  with  a quasi-symmetric  structure  that  allows  for  separation 
of  spatial  and  temporal  variations  which  generally  are  indistinguishable  in 
measurements  from  a single  satellite.  Survey  plots  for  this  event  are  shown 
in  Figure  7 in  a format  similar  to  Figure  1.  One  sees  significant  fluxes  of 
upflowing  ions  and  substantially  widened  loss  cones  in  CMF.B  in  the  central 
portion  of  the  event  (spins  2,  3,  and  4)  indicative  of  a potential  drop  be- 
low the  spacecraft.  The  maximum  energy  at  which  the  upstreaming  ion  fluxes 
were  observed  was  1.76  keV  on  spin  3.  In  this  sane  region  one  sees  90“ 
peaked  pitch  angle  distributions  with  an  up-down  symmetry  in  CMEA  similar 
to  those  found  in  the  low  energy  channels  in  the  July  29  event.  There  are 
also  90®  minimums  in  the  OlEB  distributions  which  are  most  obvious  in  the 
wings  of  the  event,  outside  the  region  of  the  90®  maxi mums  in  the  lower 
energy  electrons.  Narrowly  collimated,  field-aligned,  dovnward-f lowing  elec- 
tron "spikes"  are  seen  at  the  edges  of  the  event  in  CMEA.  litis  is  a com- 
■only  observed  signature  in  this  channel  at  the  edges  of  the  ionospheric 
structures.  We  interpret  these  data  to  indicate  that  there  is  a field- 
aligned  potential  drop  above  the  spacecraft  throughout  the  event,  with 
weak  4^  nt  the  edges  of  the  event  where  CMEA  is  responding  to  primary  elec- 
trons in  the  field-aligned  spikes;  and  an  increasing  $ toward  the  center  of 

A 

the  event  where  the  depression  in  the  CMEA  flux  levels  and  their  up-down 

symmetry  implies  that  $ > .24  keV. 

A ^ 
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The  symmetric  nature  of  the  ionospheric  structure  in  this  event 
is  best  illustrated  by  intercomparing  the  pitch  angle  distributions  at 
the  high  and  low  latitude  edges  of  the  event  by  "folding"  Figure  7 so 
that  the  two  CMEA  spikes  are  superimposed.  This  has  been  done  in  Figure 
8 where  both  the  CMEA  and  CMEB  responses  are  shown  plotted  versus  pitch 
angle  rather  than  time  in  order  to  provide  a more  precise  intercompar ison . 
As  Indicated  in  the  figure,  time  progresses  from  left  to  right  for  the 
data  from  the  high  latitude  edge  of  the  structure  and  from  right  to 
left  for  the  data  from  the  low  latitude  edge.  The  remarkable  correspon- 
dence indicates  that  we  are  traversing  a temporally  stable  structure 
and  that  the  observed  variations  are  primarily  spatial  and  angular.  It 
also  suggests  that  there  is  some  significance  to  this  pattern  at  the 
edges  of  the  structure  in  terms  of  the  mechanism  responsible  for  its 
origin. 

A CMEB  pitch  angle  distribution  measured  near  the  center  of  the 
event  is  shown  in  Figure  9.  As  in  the  previous  example,  the  90° 

■inimura  is  possibly  somewhat  obscured  by  a contribution  from  fluxes 
at  the  low  end  of  the  energy  window  with  90°  peaked  distributions  similar 
to  those  observed  in  CMEA.  Again  this  seems  plausible  since  it  is 
unlikely  that  the  transition  from  90°  maxiraums  to  90°  rainimums  in  the 
pitch  angle  distributions  occurs  exactly  at  the  edge  of  the  energy 
window.  Also,  similar  to  the  July  29  example,  the  uncertainty  in  ac 

introduced  by  such  a 90°  peaked  contribution  is  small  (see  Figure  9). 

• % 

From  the  up-down  symmetry  of  CMEA  we  can  infer  that  all  the  primaries 
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have  been  accelerated  out  of  its  passband.  This  sets  a lower  limit  on 
♦A  * .24  keV.  The  low  response  in  CMEC  suggests  that  the  spectrum  is 
falling  in  the  CMEB  energy  range,  so  for  a lower  limit  on  the  parameter 
d we  use  the  midpoint  of  the  CMEB  energy  band  for  E and  obtain  d j 1200  km 

O 

through  t he  use  of  Equation  2. 

As  bef  ore  an  alternative  limit  can  be  obtained  for  models  represented 
by  the  family  of  curves  between  the  two  extremes  shown  in  Figure  9.  For 
these  cases,  and  for  the  upper  curve,  where  there  is  a residual  response 
at  90°  due  to  the  assumed  isotropic  incident  spectrum,  one  can  set  Eg  in 
Equation  2 equal  to  the  upper  edge  of  the  CMEB  energy  window  (1.13  keV) 


and  obtain  a lower  limit  on  B<./B  j corresponding  to  d 4 1000  km. 

Burch  et  al.  [1976]  have  found  from  measurements  on  a low  altitude 
satellite  beneath  the  electron  acceleration  regions  in  "inverted  V" 
events  that  the  electron  distribution  functions  are  well  described  by 
Maxwell  ian  primary  electron  beams  which  have  been  accelerated  through 
an  electrostatic  potential.  A Maxwellian  distribution  has  the  property 
that  acceleration  through  an  electrostatic  potential  changes  the  magni- 
tude but  not  the  shape  of  the  energy  spectrum  for  energies  greater  than 
the  value  of  the  electrostatic  potential.  Thus  under  this  assumption, 
the  OlEB/CMEC  response  ratio  R is  a unique  measure  of  the  temperature  of 
the  Maxwellian,  independent  of  $ , as  long  as  $ is  less  than  the  lower 

edge  of  the  CMEB  energy  window  or  .35  keV.  As  | increases  above  this 

A 

level  it  will  depress  this  ratio  since  the  primary  electrons  arc  excluded 
from  the  lower  portion  of  the  CMEB  window.  A plot  of  the  CMEH/CMEC  response 
ratio  for  this  event  is  shown  in  Figure  10.  The  CMEB  response  was  taken 
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at  the  peak  of  the  pitch  angle  distribution  on  each  spin  while  the  CMEC 
response  was  averaged  over  the  approximately  isotropic  region  outside  the 
loss  cone.  A plausible  interpretation  of  these  data  is  that  the  tempera- 
ture is  slowly  varying  throughout  the  event  (except  for  the  region  of  ob- 
vious hardening  near  spin  5+  associated  with  an  increase  in  the  CMEC  response 
see  Figure  8)  following  the  dashed  curve  in  Figure  10.  The  depression  in  R 

from  the  dotted  curve  is  then  due  to  increasing  above  .35  keV.  This  inter- 

A 

pretation  is  supported  by  the  fact  that  we  know  from  the  depression  in  CMEA 
that  is  in  fact  increasing  toward  the  center  of  the  event,  and  by  the 
approximate  constancy  in  CMEC  which  would  most  likely  be  affected  by  a 
change  in  temperature.  Under  this  interpretation,  for  spin  3+  we  obtain 
♦A  = .54  keV  from  the  measured  depression  in  R. 

The  widest  CMEB  loss  cone  in  this  event,  indicative  of  the  maximum 
value  of  occurred  on  spin  3+  and  is  shown  in  Figure  9.  Using  this 
value  and  the  E^  obtained  from  fitting  the  above  described  Maxwellian 
to  the  CMEB  response  function,  we  obtain  $ = 1.9  keV  from  Equation  1. 

D 

Thus  the  total  parallel  potential  drop  in  this  case  is  approximately  2.4 

keV. 

Another  useful  inference  about  the  geometry  of  this  event  is  obtained 
from  the  fact  that  the  90®  minimums  in  CMEB  are  wider  at  the  edges  of  the 
event  (Figure  8)  than  in  the  center  (Figure  9).  We  infer  from  the  fact 
that  the  CMEA  response  has  not  yet  fallen  to  the  intensity  level  of  the 
secondaries  and  degraded  primaries  that  $ is  less  in  these  edge  regions  than 

A 

in  the  central  portion  of  the  event.  For  a constant  d , the  90®  minimums 
should  get  narrower  as  $A  gets  smaller.  Therefore  d must  be  less  at  the 
edges  than  in  the  center. 


A.  Narrow  spatial  regions  wun  muas 


10^  mV/m  which  arc  referred  to  as  electrostatic  shocks  by  the 

ers  . 

Berkeley  experiment/  [Mozcr  ct  al.,  1977;  Torbcrt  and  Mozer,  1978J 
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IV.  SUMMARY  AN1)  CONCLUSIONS 


We  have  shown  two  examples  of  plasma  measurements  from  within  large 
scale  ionospheric  structures  containing  parallel  electric  fields.  In  the 
July  29,  1976  example  there  is  an  "inverted  V"  electron  event  of  4°  lati- 
tudinal width  within  which  there  are  two  broad  regions  of  parallel  electric 
fields  and  a number  of  intense,  narrow  electrostatic  shocks.  Region  I is 
inferred  to  have  a parallel  potential  drop  above  the  spacecraft  of  5 kV 

3 

or  greater  with  a vertical  dimension  extending  to  > 10  km  above  the  space- 
craft. Region  II  has  a similar  field  geometry  above  the  spacecraft  and 
in  addition  shows  evidence  of  a potential  drop  below  the  spacecraft  of 
13  kV  or  greater.  In  the  September  15,  1976  example  we  infer:  1.  A 

stable  spatial  structure  of  approximately  2°  latitudinal  width;  2.  A 

3 

vertical  dimension  extending  to  > 10  km  above  the  spacecraft,  larger 

in  the  center  of  the  structure  than  at  the  edges;  3.  A total  field  aligned 

potential  drop  of  2.4  kV  with  both  i}1  and  - reaching  their  maximum  values 

in  the  center  of  the  structure.  - All  of  these  features  are  qualitatively 

consistent  with  the  V-shaped  potential  models  proposed  by  Ournett  [1972], 

Swift  [1975],  and  others  to  explain  "inverted  V"  events. 

The  simultaneous  observation  of  signatures  of  » and  $ implies  that 

the  electric  fields  involved  extend  over  a large  vertical  dimension.  If  we 

assume  that  the  fields  determined  from  the  inferred  potentials  and  vertical 

scale  sizes  are  continuous  and  roughly  independent  of  altitude,  we  obtain 

field  strengths  of  the  order  of  ■v  mV/m  in  these  broad  scale  regions  in  con- 
2 

trast  to  the<vl0  mV/ra  fields  reported  in  the  narrow  electrostatic  shocks 
(Hozer  et  al.,  1977].  The  inferred  vertical  scale  sizes  are  much  larger  than 
characteristic  dimensions  such  as  the  doByc  length  or  the  ion  gyro  radius; 


over  vertical  dimensions  of  % 103  km.  These  are  possibly  a different  class  of 
phenomena  than  those  being  studied  by  the  Berkeley  group  [Mozer  et  al. , 1977] 
whose  experiment  more  naturally  focuses  on  high  electric  field  phenomena 
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however,  this  analysis  cannot  in  fact  differentiate  between  extended  regions 
of  low  field  and  multiple  double  layers  with  regions  of  intense  fields  ex- 
tending over  short  vertical  dimensions  situated  both  above  and  below  the 
satellite. 
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FICURE  CAPTIONS 

Figure  1 - Survey  plots  for  July  19,  1976.  Universal  time  (labeled 
SYST)  is  given  on  the  abscissa.  The  location  of  four  electrostatic 
■hocks  reported  by  Mozer  el.nl.  (1977)  arc  indicated  by  broad  horizontal 
lines  along  the  abscissa.  The  energy  setting  of  the  ion  mass  spectrometers 
(step  number)  is  indicated  at  the  top  (see  Table  I.). 


Figure  2 - Expected  electron  pitch  angle  distributions  in  the  presence 
of  a parallel  electric  field. 


Figure  3 - OlED  angular  distribution  on  spin  3t  July  29,  1976.  See 
figure  1 for  definition  of  spin  numbers. 

Figure  A - CMEA,  CMEU  and  CMEC  angular  distribution  on  spin  3+, 

July  29,  1976. 

Figure  5 - CMED  angular  distribution  on  spin  7+,  July  29,  1976. 
Figure  6 - CMEC  angular  distribution  on  spin  7 and  8,  July  29,  1976. 


Figure  7 - Survey  plots  for  September  15,  1976.  Format  is  similar  to 

figure  1 . 


Figure  8 - CMEA  and  CMEB  angular  distributions  at  the  two  edges 
of  the  acceleration  region.  The  vertical  arrows  indicate  the  direction 
of  motion  of  the  measured  electrons. 


A-21 


i potential  drops  (c.g.  , 
i0“'>1C  °bU"“  Sh'°'kS-  — — resistivity)  vhich  have  chiefs- 


IMSC/D673078 


o 


Figure  9 - C1JEB  angular  distribution  on  spin  3+,  September  15,  1976. 
See  figure  8 for  definition  of  spin  numbers. 


Figure  10  - CMEB  to  CMEC  response  ratios  over  the  acceleration 
region  on  September  15,  1976.  The  right  hand  ordinate  shows  the 
temperature  of  a Maxwellian  determined  by  this  ratio. 
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TABLE  I.  DETECTOR  CHARACTERISTICS 


DETECTOR 


PARTICLE 


ENERCY,  kcV 


Energy  per  unit  charge 
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Scale  of  Electric  Field  Alone,  Magnetic 
Field  in  an  Invertcd-V  Event 
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The  characteristics  of  magnetic-f ield-aligned-electrostatic-potential 
differences  in  the  structure  of  an  invcrtcd-V  event  were  examined  using  data 
on  the  pitch-angle  and  energy  distributions  of  electrons  and  ions  obtained 
with  the  polar-orbiting  S3-3  satellite.  Potential  differences  were  found  to 
occur  simultaneously  above  and  below  the  altitude  of  the  satellite  (=7260  km) . 
Both  of  these  potential  differences  were  determined  from  an  analysis  of  the 
electron  distributions.  The  analysis  also  yielded  the  magnetic  field  inten- 
sity, B^,  above  the  satellite  where  the  potential  was  zero.  This  value  of 
B^  indicates  that  the  potential  extends  above  the  satellite  a distance 

*1200  km,  implying  an  apparent  electric  field  of  = 1 rnV/ra.  The  total  potential 

difference  was  found  to  have  an  invcrted-V  distribution  in  latitude,  with  a 
peak  value  of  =4.3  kV  and  a latitudinal  width  of  =1°.  Upward-flowing  and 

0+  ions,  with  pitch  angles  =10°  were  observed  in  the  structure.  Their  energy 

spectra  were  peaked  at  energies  near  the  potential  differences  below  the 
satellite,  but  were  very  broad  indicating  perpendicular  heating  and  an  ion 
source  within  the  potential  region.  Implications  to  theory  are  discussed. 
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Moreover,  analysis  of  the  magnetometer  data  obtained  on  the  S3-3  satellite 

during  several  traversals  of  the  potential  region  revealed  field-aligned 

2 

current  densities  in  those  regions  of  1-10  UA/m  (.Torbert  and  Jlozcr,  1977)  . 

The  observations  indicate  that  a current-driven  instability,  generating 
electrostatic  waves,  may  provide  the  mechanism  necessary  to  support  the  field- 
aligned  potential  differences  in  the  absence  of  collisions  (see  review  by 
F31thamnar , 1977) . Three  such  instabilities  have  been  proposed,  which  may  be 
identified  experimentally  by  the  scale  of  the  electric  field:  (i)  the  electric 
double  layer  (Block , 1975 ; Block  and  Fa It hammer , 197G)  has  a thickness  measured 
in  Debye  lengths  ( < 10m  at  an  altitude  of  1 in  the  auroral  zones) ; (ii)  the 
oblique  electrostatic  shock  ( Swi f t , 1975)  has  a thickness  measured  in  the  cyro- 
radii  of  the  energetic  ion  component;  and  (iii)  anomalous  resistivity  (Kir.del 
and  Kennel,  1971)  predicts  a scale  >A J/Ew»  where  A i>  is  the  potential  difference 

along  the  magnetic  field  and  E is  the  electric  field  of  the  wave  turbulence 

w 

(Falthammar , 1977) . 

In  this  report  we  will  present  an  example  of  the  ion  and  electron  data 
obtained  by  the  Lockheed  instruments  (see  e.g.  Shelley  et  al . , 1976)  on  a 
Satellite  pass  through  one  of  the  potential  regions.  The  characteristics  of  the 
potentials  inferred  from  these  measurements  will  be  described,  and  the  implicatic 
of  the  results  to  theory  will  be  discussed. 
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MEASUREMENTS 

• % 

The  Lockheed  experiment  (see  She  I ley  ct  al . . 1076;  slum  et  al.  ,‘1977) 
on  the  S3-3  satellite  measures  the  energy  and  pitch-angle  distributions  of 
ions  in  the  cnergy-pcr-charge  range  0.S  to  16  keV  and  electrons  in  the  range 
70  eV  to  24  keV.  The  measurements  are  obtained  with  three  ion  mass  spectro- 
meters and  four  magnetic  electron  spectrometers  mounted  such  that  their  view 
directions  are  perpendicular  to  the  spin  axis  of  the  satellite.  The  spin  rate 
of  the  satellite  is  ^ 3 RPM  about  an  axis  perpendicular  to  the  orbital  plane; 
hence  a nearly  complete  pitch-angle  scan  is  obtained  in  about  10  sec.  The  ion 
spectrometers  sample  the  mass-pcr-cliarge  (M/Q)  distribution  in  the  range  1-30 
once  per  second.  Each  ion  spectrometer  has  4 energy-per-charge  settings  which 
are  stopped  every  16  sec.  The  energy  settings  of  each  ion  spectrometer  to- 
gether with  the  energy  response  of  the  electron  spectrometers  are  listed  in 
Table  I. 

A survey  plot  of  the  data  obtained  over  the  northern  hemisphere  auroral 
zone  on  September  15,  1976  is  shown  in  Figure  1.  At  the  bottom  of  the  chart 
*re  shown  the  universal  time  (SYST) , longitude,  latitude,  altitude,  invariant 
latitude  (ILA),  and  magnetic  local  time.  The  four  lowest  panels  show  the 
logarithm  of  the  counts  per  half-second  of  the  electron  spectrometers.  The 
panel  labeled  PITCH  shows  the  pitch  angle  of  the  measured  particles.  Here,  the 
pitch  angle  of  0°  denotes  particles  moving  directly  down  the  field  lines.  The 
next  four  panels  show  the  logarithm  of  the  sum  of  the  counts  per  second  of  the 
3 mass  spectrometers  for  ions  of  M/Q  e 1,  2,  4,  and  16,  respectively.  A code 
designating  the  energy  steps  of  the  mass  spectrometers  is  in  the  top  panel. 
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Although  upward  moving  ions,  closely  aligned  with,  the  magnetic  field,  are 
good  indicators  of  a potential  difference  below  the  satellite,  they  do  not 
provide  accurate  information  on  the  magnitude  of  the  potential  difference; 
their  energy  spectra,  as  discussed  by  C.hieln-etti  et  al  (1078)  , are  not  simply 
related  to  the  potential  drop. 

The  satellite  entered  the  potential  region  following  the  spin  period 
denoted  by  the  number  1 under  the  pitch-angle  panel  (1052:47  UT)  and  left  the 

region  on  spin  5 (1054:37  \JT) . The  spin  numbers  are  labeled  at  the  centers 

o 

( a = 100  ) of  the  spin  periods.  Note  that  for  many  spin  periods  prior  to  1, 
the  counting  rates  of  the  electron  detectors  were  fairly  uniform  except  for 
the  loss  cones  due  to  the  atmosphere  below  the  satellite.  In  this  region,  the 

counting  rates  imply  a Maxwellian  electron  distribution  with  a number  density 

3 

of  0.66/cm  and  a temperature  of  0.8  keV.  On  spins  2 and  3 note  that  all  the 
loss  cones  become  wider  and  deeper,  those  of  the  lower  energy  electrons  being 
-■ore  strongly  affected;  the  CME  B counting  rate  outside  the  loss  cone  increases 
somewhat;  and  the  CME  A counting  rate  becomes  lower  and  has  maxima  at  pitch 
angles  near  90°.  These  characteristics,  as  discussed  below,  are  due  to  a 
potential  difference  along  the  magnetic  field,  principally  below  the  satellite, 
accelerating  electrons,  downward.  On  spin  4 the  counting  rates  of  the  CME  B 
and  CME  C detectors  develop  "butterfly'1  pitch-angle  distributions  (minima  at 
90°).  These  distributions  imply  a much  higher  potential  difference  above  the 
satellite.  Note  from  the  M/Q  = 1 and  16  panels  that  )l+  and  0+  ions  were 
moving  upward  along  the  field  during  spins  2,  3,  and  4. 
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ANALYSIS  OF  ELECTRON  MEASUREMENTS 

The  potential  differences  along  the  magnetic  field  can  be  determined 

from  the  widths  and  depths  of  the  loss  cones.  At  two  points,  (B  ,$  ) and 

s s 

(*»♦)»  along  the  dynamical  trajectory  of  an  electron,  the  energy,  w,  and 
pitch  angle,  ct,  of  the  electron  are  related  by  the  equations  expressing  the 
conservation  of  the  magnetic  moment  of  the  electron  and  the  total  energy, 

viz., 


s , 2 w 2 

sin  a = — sin  a 
D S D 

8 


w « w + e($-$  ) 
s s 


Here,  B is  the  magnetic  field  intensity,  $ is  the  electrostatic  potential, 
and  e is  the  absolute  value  of  the  electron  charge.  The  subscript,  s,  denotes 
the  values  of  the  parameters  at  the  location  of  the  satellite.  Hence,  if 
Bt  designates  the  limiting  field  intensity  at  the  "top"  of  the  atmosphere  where 

the  electron  mirrors  before  suffering  collisional  effects,  and  where  $ = S'  , 

m 

the  edge  of  the  loss  cone,  = a^,  is  given  by  the  equation 


e(<>  -»  ) 
m s 


The  potential  toward  lower  altitudes,  is  assumed  to  increase  monotonically 
from  0 to  its  maximum  value,  * , and  then  to  remain  constant.  The  affect  of 

ID 

the  shape  of  ^ along  the  magnetic  field  is  mentioned  below  in  the  discussion 
of  Figure  6. 
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Expanded  plots  of  the  counting  rates  of  the  CUE  B and  CME  C detectors 
•re  shown  in  Figure  2 and  3.  Here,  the  spin  numbers,  with  the  superscripts  - 
•nd  + denoting  respectively  the  pitch  angle  scans  obtained  during  the  first 
and  second  halves  of  the  spin  period,  are  shown  on  the  curves.  These  figures 
•how  clearly  the  increase  in  the  widths  and  depths  of  the  loss  cones  from 
•pin  1 to  spin  3.  However,  the  value  of  obtained  throuch  the 

ns  ° 

applications  of  Eq.  (2)  is  somewhat  uncertain,  principally  because  the 
detectors  respond  to  a fairly  wide  range  of  energies,  rendering  an  uncertainty 
^ in  the  effective  value  of  wg  and  smearing  ac . 
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An  improved  determination  of  ^ ~ *s  was  obtained  by  computing  the 
energy  and  pitch-angle  distributions  of  electrons  that  were  bacLscattered 
by  the  atmosphere  using  the  AURORA  code  (Walt  et  al  . . 19GR) . The  backscattered 


flux,  computed  as  discussed  bolow,  was  integrated  over  the  energy  response 
functions  of  the  detectors  to  determine  the  expected  counting  rates  of  the 
detectors  in  the  region  of  the  loss-cone  "knee"  for  various  potential  dif- 
ferences below  the  satellite.  The  value  of  f - <•  was  then  determined  by 

ra  s 

comparing  the  locations  of  the  computed  knees  with  the  experimental  values  at 
the  points  of  highest  curvature.  The  AURORA  code  was  also  used  as  described 
below  to  determine  the  potential  difference  from  the  deepending  of  the  loss 
cones.  (The  loss  cones  deepen  because  only  those  backscattered  electrons  with 


sufficient  energy  to  penetrate  the  potential  - $ can  reach  the  satellite.) 

m s 


Many  of  the  backscattcrd  electrons  are  reflected  by  the  potential  «.*  and 

m 


returned  to  the  atmosphere  where  they  combine  with  the  primary  electrons  to 
enhance  the  backscatter  (Kvans , 1975) . This  effect  was  taken  into  account  by 
using  an  iteration  procedure  to  compute  the  equilibrium  flux.  The  primary  flux 
«t  B^,  the  magnetic  field  intensity  where  < = 0,was  assumed  to  be  the  Maxwellian 

i 

distribution  that  was  measured  on  spin  1.  This  flux  was  transformed  through 

to  B^,  using  Liouvillc's  theorem,  where  it  was  used  as  input  to  the  AURORA 

code  to  compute  the  first  approximation  to  the  backscattered  flux.  Those 

•lcctrons  which  were  reflected  by  the  potential  $ were  then  added  to  the  pre- 

in 

viously-dctermined  incident  flux,  and  this  flux  was  again  used  as  input  to  the 
AURORA  code  to  compute  the  next  approximation  to  the  backscattered  flux.  The 
computation  was  continued  in  this  manner  until  the  backscattered  flux  converged 


to  • steady  value.  Tho  flux  was  then  transformed  from  B to  B , through  the 

t s 


» 
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potential  difference  4 - 4>  , to  compute  the  associated  counting  rates  of 

JO  s * 

the  detectors  in  the  loss  cone.  The  computations  were  performed  for  values 

~ *5  of  ^*5  » 3,  and  4.5  kV.  Figure  4 shows  the  corputcd  widths  of 

the  loss-cones  for  the  CME  B and  CME  C detectors,  measured  to  the  largcst- 

curvaturc  point  of  the  knee,  plotted  against  4 - 4 . These  results  are  not 
yery  sensitive  to  the  distribution  of  the  clcctYonsSnt  4 = 0 for  tn-*.s  >1.5  kV. 

In  the  center  of  the  loss  cones,  the  computed  counting  rates  were  less  than 

the  actual  counting  rates  of  the  detectors,  even  for  spin  1 where  the  potential 

was  essentially  zero.  This  discrepancy  may  exist  because  the  version  of  the 

AURORA  code  used  for  this  analysis  does  not  compute  the  production  of  secondary 

electrons,  which  may  contribute  significantly  to  the  backscattered  flux.  Because 

of  this  difficulty,  the  loss-cone  depth  computed  for  each  value  of  4 -4  was 

m s 

divided  by  the  depth  computed  for  4 -4  =0,  and  these  ratios  were  corr.nared 

m s 

with  the  corresponding  experimental  ratios,  i.e.,  the  loss-cone  depth  observed 

in  the  potential  region  divided  by  the  depth  observed  on  spin  1.  Here,  the 

loss-cone  depth  is  defined  as  the  ratio  of  the  counting  rate  at  to  the  counting 

rate  at  a = 100  . This  procedure  tends  to  cancel  errors  in  the  backscattered 

flux  which  are  proportional  to  the  flux.  The  computed  ratios  of  the  loss-conc 

depths  for  the  CME  C detector  are  shown  in  Figure  5 plotted  against  4 - 4 . 

m s 

The  resulting  potential  differences  below  the  satellite,  inferred  from 

both  the  widths  and  depths  of  the  loss  cones  on  the  pass  through  the  potential 

region,  arc  shown  in  Figure  6.  The  values  determined  from  the  widths  of  the 

loss  cones  arc  regarded  to  be  the  more  accurate,  and  the  uncertainties  of 

these  values  due  to  the  unccrtaint ics  in  a arc  estimated  to  be  within  the 

c 

tttoT  bars  shown  in  the  figure.  Another  source  of  error  may  be  due  to  an 
effect  of  the  shape  of  the  potential  distribution  In  the  calculation  of  the 
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electron  trajectory  will  go  through  those  points.  The  electron  moving 
upward  from  the  atmosphere  could  still  be  reflected  backward  before  reaching 
Bs  because  of  an  appropriately  stiaped  potential  distribution  between  the 
points  (sec  Chiu  and  Schulz,  1978;  Whipple . 1977;  and  Cladis,  et  al . . 1977). 

Chiu  and  Schulz  (1978)  point  out  that  such  reflection  docs  not  occur  if  the 
potential,  ♦(&),  when  plotted  against  B,  is  concave  downward  between  the 
two  points.  In  general,  the  effect  is  less  likely  to  occur  if  the  scale  of 
the  electric  field  is  large;  according  to  Cladis.  ct  al.  (1977),  a necessary 
condition  for  the  reflection  is  that  the  downward  electric  force  eV^i  acting 
on  the  electron  be  greater  than  the  oppositely-directed  magnetic  force 
at  least  at  one  point  within  the  interval  (p  is  the  magnetic  moment  of  the 
electron,  and  V is  the  component  of  the  gradient  parallel  to  the  magnetic 
field).  This  is  not  a sufficient  condition,  however:  Even  if  £(B)  is  a step 
function,  it  can  still  be  penetrated  by  an  electron  that  satisfies  (1)  if  the 
potential  step  is  located  at  a sufficiently  low  value  of  B within  the  interval. 
Nevertheless , two  aspects  of  the  data  indicate  that  an  appreciable  error  was 
not  introduced  by  the  potential-shape  effect:  (i)  for  each  spin  period,  the 
potential  differences,  ♦ , inferred  from  the  edges  of  the  pitch-angle  distri- 

butions given  by  the  counting  rates  of  the  detectors  with  different  pass-band 
energies  (CME  B and  CME  C) , were  in  good  agreement;  and  (ii)  the  shapes  of 
the  counting-rate  curves  for  the  different  detectors  in  the  regions  of  the 
pitch-angle  knees  were  quite  similar  to  the  corresponding  counting-rate  curves 
based  on  the  computed  backscattcrcd  flux. 
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The  potential  differences  $ above  the  satellite  are  also  shown  in 

6 

Figure  6.  The  value  of  $ on  spins  2 and  3 was  estimated  from  the  increase 

S 

in  the  counting  rates  of  the  CME  B and  CME  C detettors  outside  the  loss 
cones,  and  the  reduction  of  the  CME  A counting  rate.  The  potential  $ 

s 

shifts  the  spectrum  and,  according  to  Liouvillc's  theorem,  increases 

the  d Erectional  flux  observed  by  the  detectors  (w  > e i ) by  the  factor 

s s 

(w  + e $ )/w  . The  value  51  *v  0.25  kV  was  estimated  from  the  increase  in 
s s s s 

the  electron  flux  implied  by  the  counting  rates  of  the  detectors  relative  to 

those  on  spin  1 and  using  for  w^  the  central  band  pass  energies.  By  falling 

through  this  potential,  the  primary  electrons  acquired  energies  beyond  the  pass 

band  of  the  CME  A detector.  The  magnitude  of  its  counting  rate  and  the  approx- 

o 

imatc  symmetry  of  the  counting  rate  about  = 90  indicate  that  it  was 
detecting  principally  backscattcred  electrons  that  were  reflected  downward  by 
the  electric  field  above  the  satellite. 
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On  the  pitch-angle  scans  2 and  4 , the  CME  A detector  was  probably  also 
counting  electrons  tliat  were  trapped,  at  least  temporarily,  between  the  mag- 
netic ficl^l  below  the  satellite  and  the  electric  field  above.  Electrons  may 
become  locally  trapped  in  this  manner  by  time  and/or  space  variations  of  the 
electric  field  along  the  electron  trajectories.  Our  analysis  of  the  CME  A 
counting  rate,  as  well  as  the  CME  li  counting  rate  on  spin  4 which  is  discussed 
below,  reveals  that  the  local  trapping  consists  mainly  of  degraded  primaries 
and  scattered  electrons  which  had  pitch  angles  near  00°  at  the  "top"  of  the 


atmosphere,  where  B = Bfc  . The  trapping  may  occur  for  energies  in  the  range, 


e (*m  - V 

Bt  ,•  2 , 

— Sin  a - 1 

B s 


< w < 
— s — 
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B1  2 
1 - Sin  a 
B s 
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This  range  increases  from  zero  at 


* « Sin'1  r Bs  ♦„ 1 ' 

[Bj  (*  - -s>  ♦ Bt  J 


(3) 


(4) 


where  the  limiting  energy  w^  ; w^  is  obtained  by  putting  = ng 
in  the  equation  on  either  side  of  (3),  to  a maximum  at  a = 90°. 
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The  right* hand  side  of  (3)  follows  from  the  conservation  of  the 


magnetic  moment  for  electrons  mirroring  above  the  satellite,  at  B where 

1 


^ = 0.  The  left-hand  side  of  (3)  follows  from  the  conservation  of  the  mag- 
netic moment  for  electrons  mirroring  at  B^.  The  minimum  energy  for  trapping 
is  given  by  the  left-hand  side  of  (3)  for  ag  = 90°.  This  energy  is  within 
the  upper  portion  of  the  pass  band  of  the  CME  A detector.  The  trapping  of 


these  electrons  may  therefore  account  for  the  sharp  peaks  of  the  CME  A counting 


rate  centered  at  - 90  on  scans  2 and  4 . Unfortunately,  the  maximum 


energy  of  the  trapped  electrons,  given  by  the  right-hand  side  of  (3),  is  much 


beyond  the  pass  band  of  the  detector  for  all  values  of  a > a , hence  this 

s s 


detector  cannot  provide  information  on  B^. 


The  CME  B and  CME  C counting  rates  on  spin  4 provided  information  on  both 
♦ and  B . Since  the  potential  above  the  satellite  increases  the  energy  ccmoor.er 
of  the  primary  electrons  along  the  magnetic  field,  the  electrons  become  more 
aligned  with  the  magnetic  field, and  the  alignment  is  closer  for  the  lower-energy 
primaries.  Such  downgoing  primaries,  together  with  those  which  are  reflected 
by  the  maqnetic  field  below  the  satellite,  form  a butterfly  pitch-angle 
distribution.  The  butterfly  shape  of  the  CME  C counting  rate  (see  Figure  3) 
on  spin  4 is  due  mainly  to  the  acceleration  of  the  primaries.  A potential  $ 
that  accounts  for  the  pitch-angle  distribution  of  the  primaries  detected  by 
the  CME  C detector  causes  the  lower-energy  primaries  observed  by  the  CME  B 
detector  to  appear  only  at  very  small  angles  to  the  field.  Hence,  the  energy 


et  must  have  been  within  the  upper  portion  of  the  CME  B pass  band  to  account 


for  both  the  higher  counting  rate  near  a “ 0 (due  to  downgoing  primaries)  , 
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relative  to  that  near  ag  «*  180°  (backscattered  electrons),  and  the  symmetry 

of  the  counting  rate  about  90°  at  pitch  angles  that  are  more  inclined  with 

the  magnetic  field.  The  latter  is  due  to  backscattered  electrons  that  were 

also  reflected  downward  by  $s»  and  locally-trapped  electrons. 

The  values  of  4g  and  B,  were  estimated  from  the  CME  C counting  rate  by 

assuming  the  primary  electrons  at.  B^to  be  isotropic  in  the  downward  hemisphere 

and  to  have  the  Maxwellian  distribution  that  was  measured  on  spin  1.  The 

directional  flux  of  these  electrons  at  the  satellite,  B , from  Liouville's 

s 


theorem,  is 


j <V  “s'  V = 


n w 

o s 


/■>  3 3.H 

(2  ir  m w ) 
e 


Exp 


[- 


(w  - c 4 )/w 
s s e 


1 


(5) 


where,  from  the  conservation  of  the  magnetic  moment. 


a **  Sin 
s 


--  [ 


(1  - e !s  ) Bs  2 

V 57  Sln 

S 1 


(6) 


for  w >_  e 4 . As  mentioned  previously,  n ^ 0.7/cm  and  w £ 0.8  kcV; 

S __l  S 0 6 

Qj  is  the  pitch  angle  of  the  electrons  at  . The  counting  rate  of  the  CME  C 
detector  due  to  this  flux  is  given  by  the  integral. 


CR  (ag)  =•  G 


/ 


5 keV 

3 <V  “s'  V dws 


(7) 


•6  2 

G ■ 5.4  x 10  cm  sr  is  the  geometric  factor  of  the  detector,  and  w^  is  the 

larger  of  W = 1.6  keV  (the  low-energy  limit  of  the  detector  band  pass)  or 

1* 


W1 


e 4 


1 - Sin2  a 

B. 


(8) 
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which  follows  from  (6)  for  = 90°.  Note  that  the  flux  given  by  (5)  is 
constant  over  the  pass  band  of  the  detector  in  the  pitch-angle  range  0 < 

.5.  where  follows  from  (8)  for  w^  a w^.  Hence,  the  counting  rate 

of  the  detector  should  be  constant  in  that  pitch  angle  Interval.  The  situa- 
is  illustrated  in  Figure  7,  which  depicts  the  response  of  the  detector, 
as  a function  of  pitch  angle,  on  the  contours  of  constant  flux  obtained  by 
transforming  the  isotropic  flux  from  B^ , 4 = 0 to  Bg , 4g.  The  relationship 
between  and  B^  at  the  limit  of  the  constant  counting-rate  interval, 
given  by  the  equation, 


* #s  = WL  (1  “ — Si°2  aL>  (9) 

B 

s 

is  valid  for  an  isotropic  flux  at  B^ , regardless  of  the  energy  distribution. 

- o 

From  Figure  3,  on  spin  4 , it  appears  that  a n,  60  . At  the  higher  pitch 

L 

angles,  the  counting  rate  decreases  because  w^  increases  and  the  integration  (7) 

is  over  the  higher-energy  portion  of  the  energy  spectrum.  The  distribution  shewn 

on  spin  4 can  be  closely  matched  with  4 = 0.8  kV  and  B./B  = .77. 
r s I s 

In  order  to  test  the  sensitivity  of  these  results  to  the  assumed  conditions, 

these  calculations  were  also  done  for  Maxwellians  of  different  temperatures  and 

and  for  lower  counting  rate  minima  at  90°.  (The  counting  rate  due  to  the  primaries 

would  be  lower  there  if  the  detector  were  also  counting  backseat tered  and 

' trapped  electrons).  A distribution  similar  to  that  on  spin  4 , but  with  zero 

counting  rate  at  90°,  was  computed  for  4g  = 0.66  kV  and  B^/Bg  = .06.  Using  4g 

■ .8kcV  and  Bj/Bg  " -77,  but  changing  the  temperature  of  the  Maxwellian  from  .8 

to  .6kcV  principally  decreased  the  counting  rate  at  90°  by  about  50\;  increasing 

tho  temperature  from  .8  tolkeV,  increased  the  90°  counting  rate  by  about  43%. 

. A -variation  of  the  90°  counting  rate  of  ♦_  50%  was  found  to  correspond  to  a variation 

in  ♦ of  +5%,  and  in  B./B  of  + 5%. 

■ IS  — 
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Noto  that  the  value  «>g  £ .8  kV,  which  is  near  the  upper  limit  of  the  pass 
band  of  the  CME  B detector,  is  also  consistent  with  the  discussion  above  regard- 
ing the  counting  rate  of  that  detector  on  spin  4.  Moreover,  by  substituting  the 
values  = .8  JCV,  = 2.6  kV,  B1/Bj;  = 0.77,  and  Bs/flt  % 1/8  in  Equation  (3), 

the  energy  range  of  the  temporarily-trapped  electrons  near  the  edge  of  the  loss 
cone,  where  the  CKE  B counting  rate  is  a maximum,  is  found  to  be  0.6  to  1.3  keV. 

The  minimum  energy  of  the  downgoing  primaries  at  B is  .8  keV,  for  $ = .8  kV; 

S s 

hence,  a high  flux  of  degraded  primaries  would  return  to  the  satellite  with 
energies  in. the  trapping  range.  The  flux  at  the  maximum,  therefore,  appears  to 

consist  of  degraded  primaries  that  are  backscattcred  and  temporarily  trapped,  as 
discussed  before,  making  at  least  one  bounce  from  the  electric  field  above  the 
satellite.  The  decrease  of  the  counting  rate  toward  a = 90°  appears  to  bo  due 
to  the  decreasing  fluxes  of  successively  lower  energy  electrons  of  the  atmospheric 
back scatter  that  are  temporarily  trapped. 

As  the  boundary  of  the  potential  region  was  approached,  near  the  end  of 
the  4 scan,  ♦ decreased  to  zero,  as  evidenced  by  the  sharp  increase  in  the 
counting  rate  of  the  CM£  A detector  (note  counting  rate  "spike"  in  Figure  1) . 
Simultaneously,  as  shown  in  Figure  1,  the  flux  of  the  high-energy  primaries  seen 
by  the  CM.E  C detector  increased  sharply  to  a maximum  (see  also  Figure  3) . Sub- 
sequently, the  CME  C and  C.ME  B counting  rates  decreased  steadily,  reaching  the 
ainimum  values  at  the  end  of  the  5 scan.  There,  the  potential  difference  below 
the  satellite  also  must  have  been  near  zero  because  the  upward-moving  ions  were 
not  observed. 
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ION  MEASl'KEMENTS 


The  upward-moving  ions  observed  on  spins  2,  3,  and  4 had  pitch  angles 

less  than  about  10°.  They  consisted  of  H+  and  0+,  but  the  0+  flux  was 

lower  than  the  H+  flux  by  about  an  order  of  magnitude.  The  energy  spectra 

of  the  11  ions  on  the  3 spins,  which  include  the  data  from  the  Aerospace 

electrostatic  analyzer  (Mizera  and  Fennel.  1978)  as  well  as  the  Lockheed 

■ass  spectrometer,  are  shown  in  the  3 panels  of  Figure  8.  The  Aerospace  data 

have  not  been  corrected  for  the  small  contribution  from  0+.  The  times  of  the 

measurements  arc  shown  in  the  panels.  Because  of  the  angular  displacement  of 

the  spectrometer  apertures  about  the  spin  axis,  the  Lockheed  observations  of 

i<* 

the  source  cones  followed  by"'  Aerospace  observations  by  =5  seconds.  The  dif- 
ferences between  the  two  data  sets  seen  in  Figure  8 could  result  from  spatial 
and/or  temporal  fluctuations  during  this  interval.  The  calibrations  of  the  two 
spectrometers  were  reconciled  by  inter-comparisons  in  the  slowly  varying,  nearly 


isotropic  fluxes  of  the  radiation  belts,  but  no  attempt  has  been  made  to  unfold 
the  angular  distributions  of  the  upflowing  ions  from  the  instrument  response 
functions.  The  values  given,  therefore,  represent  averages  over  the  fields  of 
view  of  the  respective  spectrometers  (6°  x 5°  full  width  for  the  Lockheed  instru- 
nent  and  10  x 25  full  width  for  the  Aerospace  instrument).  Since  the  field  of 
view  of  the  Aerospace  instrument  was  larger  than  the  angular  width  of  the  ion 
beams,  the  Aerospace  points  should  be  considered  strictly  as  lower  limits.  The 

i 

spectral  shapes  shown  should  be  valid,  however,  if  the  cone  widths  are  not  a 
Strong  function  of  energy. 

Since  the  Lockheed  mass  spectrometer  is  designed  to  scan  rapidly  with 
respect  to  mass,  it  acquires  only  a three-point  energy  distribution  on  each  spin. 
The  more  complete  spectrums  determined  by  the  Aerospace  spectrometers  Indicate 
that  the  peak  fluxes  arc  at  energies  that  are  roughly  comparable  to  the  values  of 
• shown  in  Figure  6 at  the  appropriate  times;  i.e.,  at  points  on  the 
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abscissa  about  one-fourth  of  a spin  period  prior  to  the  numbered  spins.  The 
fewer  Lockheed  data  points  indicate  that  the  peaks  may  have  been  at  higher  energies 
on  spins  2 and  3,  and  at  a lower  energy’ on  spin  4,  in  agreement  with  Figure  6. 

Note  also  that  the  spectra  are  very  broad,  much  more  so  than  expected  for  a 
pure  acceleration  by  the  potential  difference,  notwithstanding  plausible  time 
variations  of  the  potential  and  the  long  transit  time  of  the  ions.  Moreover, 
according  to  Ghiclmetti  ct  al.  (1978),  ions  with  conical  pitch-angle  distributions, 
with  maxima  ranging  from  very  small  angles  to  90°,  are  also  occasionally  observed 
in  the  potential  regions.  It  appears,  therefore,  that  the  ions  are  being  accel- 
erated perpendicular  to  the  magnetic  field,  as  well  as  in  the  parallel  direction, 
in  the  potential  region.  Since  the  maxima  of  the  spectra  shown  in  Figure  8 are 
correlated  fairly  well  with  the  potential  differences  below  the  satellite,  a high 
fraction  of  the  ions  must  have  fallen  through  the  entire  potential  drop. 


DISCUSSION 


The  inverted-V  distribution  of  the  potential  difference  as  a function  of  spin 
number,  shown  in  Figure  6,  is  in  agreement  with  the  similarly-shaped  energy 
versus  time  displays  of  the  precipitating  electrons  (the  invcrtcd-V  electron  events 
Observed  at  lower  altitudes  (Frank  and  Ackersoii,  1972)  and  of  the  upward-moving 
ions  observed  on  the  S3-3  satellite  (Ghiclmetti  et  al.  , 1978 Minora  and  Fennell. 
1977).  As  discussed  by  Ghielmetti  et  al . , (1978),  the  potential  structure,  as 
inferred  from  the  ion  measurement,  is  evidently  a common ly-oc curing  feature  of  the 
auroral  zone.  Analyses  of  different  sets  of  the  S3-3  data,  have  yielded  results 
simil  ar  to  those  shown  in  Figure  6 (see  e.g..  Sharp  et  al.,  1978;  Mizcra  and 
Fennell , 1977) , including  simultaneous  potential  differences  above  and  below  the 
satellite . 

The  apparent  electric  field,  E^,  along  the  magnetic  field  above  the 

satellite  can  be  estimated  from  the  values  of  B,  and  ♦ obtained  from  the 

1 6 
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spin  4 data.  Tho  ratio  B /B  = .77  corresponds  to  a distance  of  about  1200 

X s 

km  above  the  satellite.  Hence  E^  £ 0.8/1200  or  .7  mV/ro.  The  value  of  E^ 
below  the  satellite  can  be  estimated  from  the  values  of  4^  - 4 by  assuming 
that  is  located  at  some  low  altitude,  say  3000  km.  Such  an  altitude  is 
reasonable  since  many  observations  of  highly  field-aligned  ( < 20°) , low- 
energy  electrons  (sometimes  fairly  monoenerget ic)  have  been  made  in  the  auroral 
cones  which  indicate  that  $ may  be  located  at  altitudes  _<  2000  km  (see,  e.g. 
Arnoldy  et  al . , 1974) ; and  Ghirlnetti  et  al . , (1978)  found  that  the  probability 
of  seeing  upward  flowing  ions  above  4000  km  cn  satellite  passes  over  the  auroral 
zones  is  high,  about  60%.  That  altitude  for  4^  implies  an  extent  of  the  potential 
below  the  satellite  of  'v  4000  km.  The  corresponding  electric  fields,  from  the 
spin  2,  3 and  4 results,  are  .4,  .9,  and  .4  mV/m,  respectively. 

The  apparent  scale  of  the  potentials  is  therefore  much  larger  than  the 
Debye  length  or  the  energetic  ion  gyroradius.  Moreover,  since  electrostatic 
ion  cyclotron  (EIC)  turbulence  is  occasionally  observed  by  tho  S3 -3  satellite 
(Mozcr  et  al  . 1977 ; Kantnei  et  al  . , 1978)  , it  ap>pears  that  the  electric  field  is  r.:;- 
likely  supported  by  anomalous  resistivity,  as  discussed  by  Kindcl  and  Kennel  (1971;. 
Papadopoulos  (1977),  and  I.ysak  et  al  . (1978). 

Such  large-scale  potential  distributions  can,  in  principle,  be  computed  from 
the  quasi-neutrality  principle  (see,  e.g.,  Pcrsson,  1963) . That  is,  the  potential 
distribution  must  be  such  that  the  electron  and  ion  number  densities  are  essen- 
tially identical  at  each  point  along  the  magnetic  field.  If  it  is  assumed  that 
tho  clcctrostatic-wavc  turbulence  prevents  the  "thermal"  ions  and  electrons 
from  moving  along  the  magnetic  field,  as  is  necessary  to  achieve  the  anomalous 
resistance,  the  charge  neutrality  must  be  established  by  the  more  energetic 
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paxticlcs.  An  analysis  by  Cl  ad  is  c t a 1 . (1978)  has  been  conducted  to 
determine  whether  the  number  densities  of  the  upgoing  ions  and  electrons 
observed  on  spin  4,  and  the  downgoing  electrons,  assumed  to  be  Maxwellian  at 

as  discussed  previously,  might  yield  such  a self-consistent  potential  distri- 
bution in  the  region  above  the  satellite  between  B,  and  B . The  fluxes  above 

1 s 

the  satellite  were  computed  from  an  application  of  Liouville's  theorem,  and 
the  potentials  at  the  boundaries  B^  and  B^  were  kept  at  the  values,  0 and  .8  kV 
respectively,  determined  from  the  electron  pitch-angle  distributions.  The 
Aerospace  data  on  the  proton  spectrum  on  spin  4 were  not  available  at  that 
time,  hence  the  analysis  was  performed  for  various  mathematical  forms  for  the 
proton  flux  that  seemed  to  bound  the  Lockheed  data  shown  in  Figure  8.  The 
proton  flux  was  normalized  such  that  the  number  density  of  the  protons  was 

equal  to  that  of  the  electrons  at  B . With  the  proton  flux  proportional  to 

-v/w  s 

P 

w e , where  w^  = .2  keV,  for  example,  excellent  agreement  of  the  number 

densities  throughout  the  interval  was  achieved  for  a potential  that  increased 

almost  linearly  with  B from  B = .785  B to  B . 

s s 

The  results  or  that  analysis  reveal  that  a self-consistent  solution  for 
the  potential  distribution  cannot  be  obtained  with  the  Aerospace  flux  shown  in 
the  third  panel  of  Figure  8.  In  the  first  place  the  flux  is  too  low  by  a factor 

i 

of  about  8 to  match  the  number  density  of  the  electrons  at  Bg.  Even  if  this 
flux  were  raised  to  match  the  number  densities  at  Bg,  the  proton  densities  would 
exceed  the  electron  densities,  regardless  of  the  potential  shape,  at  the  lower 
B values.  A self-consistent  solution  can  be  obtained  only  if,  toward  higher 
altitudes,  the  ion  spectrum  becomes  progressively  softer  and/or  the  ion  pitch- 
angle  distribution  becomes  progressively  wider.  Both  such  effects  would  be 
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expected  from  a streaming  instability  of  the  ions. 

If  the  fluxes  above  the  satellite  are  correctly  given  by  Liouville's 
theorem,  it  would  then  appear  that  lower-energy  electrons  and/or  ions,  not  in- 
cluded in  the  observations , are  contributing  to  the  necessary  transport  of 
charge.  From  the  standpoint  of  anomalous  resistivity  theory,  these  particles 
might  be  due  to  run-away  electrons  from  the  locally-heated  electron  distribution 
and/or  ions  that  are  accelerated  locally  by  resonating  with  the  electrostatic 
ion  cyclotron  waves..  In  fact,  such  ions,  which  escape  from  the  local  turbulence, 
and  are  subsequently  accelerated  by  a portion  of  the  potential  difference  below 
the  satellite,  may  also  account  for  the  ion  fluxes  shown  in  Figure  8 at  energies 
below  the  peaks.  Of  course,  .a  streaming  instability  of  the  ion  beam  might  also 
explain  the  presence  of  the  low  energy  ions. 

Quite  often,  within  the  potential  region,  one  or  more  narrow  regions 

o 

( < .1  in  latitude) appear  which  arc  characterized  by  high  electric  fields 

2 

( •u  10  mV/m),  principally  perpendicular  to  the  magnetic  field,  with  associated 
potentials  in  the  kV  range  (Mozcr  ct  al.,  1977).  These  structures  may  be  the 
oblique  electrostatic  shocks  discussed  by  Swi f t (1975) . An  observation  of 
several  such  "shocks"  within  the  potential  region  is  discussed  by  Sharp  ct  al . 
(1978).  They  might  represent  different  dynamical  states  of  the  same  current 

l 

driven  instability. 

I 
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CONCLUSIONS 


The  simultaneous  electrostatic  potential  differences  above  and 
below  the  satellite  inferred  from  the  ion  and  electron  measurements,  indicate 

I 

that  the  scale  of  the  potential  along  the  magnetic  field  is  large.  This 
characteristic,  together  with  observations  of  electrostatic  turbulence  and 
high  field-aligned  currents  by  the  S3-3  satellite,  indicate  that  the  electric 
field  along  the  magnetic  field  may  be  supported  by  anomalous  resistivity. 

A self-consistent  solution  for  the  potential  distribution  above  the 
satellite  cannot  be  found  for  the  observed  distributions  of  the  upgoing  ions 
and  electrons  and  the  assumed  J-'-axwel lian  electron  distribution  above  the 
satellite  where  the  potential  is  zero.  The  quasi-neutrality  condition  can  be 
satisfied  only  if  (i)  the  ion  distribution  is  altered  above  the  satellite  by 
interactions  that  continually  soften  the  spectrum  and/or  widen  the  pitch-angle 
distribution  or  (ii)  particles  with  energies  lower  than  those  seen  by  the 
detectors  arc  also  transported  along  the  field.  A streaming  instability  of  the 
ion  beam  might  account  for  (i) . Likely  candidates  for  (ii)  arc  particles  which 
may  escape  from  the  local  turbulence:  run-away  electrons  from  the  heated  electro 
distribution  and  ions  which  are  perpendicularly  heated  by  EIC  waves.  Such  a 
source  of  the  ions  within  the  potential  structure  below  the  satellite  is  con- 
sistent with  the  observed  energy  spectra  of  the  ions,  which  reveal  the  presence 
of  ions  with  energies  less  than  the  potential  differences  below  the  satellite. 
Nevertheless,  since  the  peaks  of  the  ion  spectra  are  at  energies  corresponding 
roughly  to  the  full  potential  differences  below  the  satellite,  it  appears  that 
most  of  the  ions  originated  near  the  altitude  of  the  maximum  value  of  the  po- 
tential. The  low-energy  ions  might  also  result  from  an  ion-streaming  instability 
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TABLE  I. 

DETECTOR 

CHARACTERISTICS 

DETECTOR 

PARTICLE 

ENERGY, 

, kcV 

2 

CDE,  cm 

sr  kc-V 

CMEA 

Electrons 

0.07  - 

0. 24 

1.2 

X 

IQ'6 

CHEB 

Electrons 

0.35  - 

1.1 

6.5 

X 

10-6 

CMEC 

Electrons 

1.6  - 

5.0 

1.9 

X 

io~5 

CMED 

Electrons 

7.3  - 

24 

6.5 

X 

io"5 

Energy  per  unit 

charge 

Step  1 

2 

3 

4 

CXA  1 

Ions 

0.50 

0.68 

0.94 

1.28 

CXA  2 

Ions 

1.76 

2.4 

3.3 

4.5 

CXA  3 

Ions 

6.2 

8.5 

11.6 

16.0 
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figure  captions 

1.  Survey  plots  for  September  15,  1976.  (see  text). 

2.  Counting  rate  of  CNE  B detector  as  function  of  pitch  angle.  The  satellite 
•pin  number,  with  the  subscripts  - and  + denoting  the  first  and  second 
pitch-angle  scans  made  during  the  satellite  spin,  arc  shown  at  the  curves. 

3.  Counting  rate  of  CME  C detector  as  function  of  pitch  angle.  (see  caption 
of  Figure  2). 

4.  Loss-cone  width  measured  to  point  of  largest  curvature  of  "knee",  as  computed 
with  the  AURORA  code,  as  function  of  potential  difference  below  satellite. 

I 

5.  Loss-cone  depths  computed  with  AURORA  code  for  various  potential  differences 
below  satellite  divided  by  the  computed  loss-cone  depth  for  the  potential 
difference  equal  to  zero,  (see  text). 

6.  Electrostatic  potential  differences  inferred  from  measurements  of  electron 
distributions  on  satellite  pass  through  potential  region.  The  circles 
and  the  triangles  denote  the  potential  differences  below  the  satellite 
determined  from  the  loss-cone  widths  and  depths,  respectively,  with  the 
error  bars  indicating  the  uncertainties  based  on  the  loss-cone  width 
determinations.  The  squares  denote  the  total  potential  differences,  which 
include  the  potential  differences  above  the  satellite. 

• \ 

Illustration  of 'detector  response  to  primary  electrons  in  downward  hemisphere. 

In  diagram  (a),  contours  of  constant  flux  at  B^,*  ■ 0 - where  the  flux  is  assumed 
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Isotropic  - arc  shown  in  velocity  space;  the  corresponding  contours  at 

B#,  ♦ are  shown  in  (b).  The  contour  segments  between  the  dots  labeled 

1,  1;  2,  2;  etc.,  in  (a)  map  into  the  similarly  labeled  segments  in  (b). 

2 

At  the  contours  are  given  the  fluxes  in  units  of  (cm  •sec-sr.kcV)  with  the 

3 

power  of  10  in  parenthesis,  for  ne  « .7/cra  and  w^  = .8  kcV  (see  Eq.  (5)), 

and  the  electron  energies  in  keV.  The  detector  response  (in  the  range  1.6 

to  5 kcV)  at  ag  is  depicted  by  the  cross-hatched  area.  Note  that  the  detector 

counting  rate  is  constant  from  a = 0 to  a = a.  . 

s s L 

Energy  spectra  of  ions  measured  by  the  Lockheed  and  Aerospace  (Mir.era  and 
Fennell , 1978)  groups  on  spins  2,  3,  and  4.  The  times  of  tiro  measurements  are 
shown  in  the  panels. 
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POTENTIAL  DITFERENCE  ALONG  MAGNETIC  FIELD 
THROUGH  STRUCTURE  OF  INVERTED-V  EVENT 


LEGEND:  £ POTENTIAL  BELOW  SATELLITE  FROM  DATA  ON  LOSS 
T CONE  WIDTHS  WITH  ESTIMATED  UNCERTAINTIES 

A POTENTIAL  BELOW  SATELLITE  FROM  DATA  ON  LOSS 
CONE  DEPTHS 

□ TOTAL  POTENTIAL  (INCLUDES  ESTIMATES  Or 
POTENTIAL  ABOVE  SATELLITE) 
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ABSTRACT 


Energetic  (keV)  ions  with  field  aligned  pitch  angle  distributions  have 
been  observed  streaming  downward  in  the  high  altitude  ('1  K ) auroral  iono- 

u 

sphere  by  the  ion-mass  spectrometer  on  the  S3-3  satellite.  These  downward 
streaming  ion  events  were  observed  much  less  frequently  titan  the  previously 

reported  upward  flowing  ion  events.  They  exhibited  peak  fluxes  of  up  to 

7 2 -1 

“7  x 10  (cm  -scc-ster-kcV)  and  typically  consisted  only  of  protons  with 
energies  less  than  '3  keV.  In  two  examples  investigated  in  detail,  the 
similarities  between  the  downward  flowing  and  coincident  trapped  ioit  distri- 
butions, suggest  a common  origin.  Spatially  localized  regions  of  enhanced 
hot  (-keV)  plasma  density  were  frequently  observed  in  the  low  L portion  of  the 
plasma  sheet.  The  statistical  location  of  these  "plasma  clouds"  correlates 
well  with  the  substorm  injection  boundary  near  dusk  inferred  by  Mcllwain  (197s) 
and  the  downward  flowing  ion  events  occur  preferentially  within  such  "plasma 
clouds".  It  is  suggested  that  upward  flowing  ions  from  the  auroral  accelera- 
tion regions  are  responsible  for  both  the  downward  flowing  ions  and  for  at 
least  some  of  the  "plasma  clouds".  It  is  concluded  that  injection  of  energetic 
ions  from  this  ionospheric  source  into  the  trapped  population  of  the  plornsa 
sheet  is  common  and  may  contribute  significantly  to  the  equatorial  plasma 
density.  The  significant  difference  in  the  frequency  of  upstreaming  and  down- 
streaming ion  events  further  suggests  that  parallel  electric  fields  involving 
potential  drops  of  J 500  volts  are  directed  preferentially  upward  in  the  altitude 


range  from  2000  km  to  ' 3 Rg. 
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A general  discussion  of  the  energetic  ion  composition  in  the  magneto- 
sphere is  provided  in  the  reviews  by  Johnson  et  al.  (1975)  and  Shelley  (1978). 
Specifically,  the  presence  of  large  fluxes  of  energetic  trapped  0+  ions  in  the 
inner  ring  current  (Johnson  et  al . , 1977)  and  of  precipitating  0+  ions  in  the 
auroral  zone  (Shelley  et  al . , 1972,  1974;  Sharp  et  al.,  1974,  1976a, b;  Johnson 
et  al . , 1975)  during  geomagnetic  storms  implies  that  a sizeable  fraction  of  the 
energetic  magnetospheric  particle  population  is  of  ionospheric  origin.  On  the 
basis  of  morphological  studies  on  these  precipitating  ions  Sharp  et  al.  (1976b) 
have  inferred  a local  injection  process  that  may  be  operating  over  a wide  region 
of  local  times.  More  recently  Shelley  et  a I . (1976)  and  Sharp  et  al.  (1977) 
have  directly  observed  energetic  ions  being  accelerated  out  of  the  ionosphere 
and  injected  into  the  magnetosphere.  Synoptic  studies  have  identified  the 
statistical  auroral  oval  and  the  altitude  range  from  2000  to  >8000  km  as  a 
source  region  where  this  upward  acceleration  of  ionospheric  ions  is  commonly 
occurring  (Ohio  1 met ti  et  al.,  1978a).  How  and  to  what  extent  this  ionospheric 
source  contributes  to  the  various  magnetospheric  particle  populations  and  its 
relation  to  other  dynamical  processes  in  the  magnetosphere  remain  open 
questions  at  present  however. 

Strong  perpendicular  electric  fields  are  frequently  encountered  in  the 
regions  of  upward  flowing  ions  (Mozer  et  al.,  1977;  Torhcrt  and  Mu/e r , 1978). 
These  field  measurements  and  the  signatures  observed  in  the  particle  distribu- 
tions have  been  interpreted  as  evidence  for  upward  directed  parallel  electric 
fields  (Mozer  et  al. , 1977;  Mizera  and  Fennell,  1977;  Torbort  and  Mozer,  1978; 
Cladls  and  Sharp,  1978;  Sharp  et  al . , 1978a).  If  oppositely  directed  parallel 
electric  fields  existed  in  these  regions  they  would  be  expected  to.  give  rise 
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to  similar  field  aligned  ion  distributions  flowing  in  the  downward  direction. 
Thus,  the  relative  frequency  of  occurrence  of  downward  flowing  ion  (Dll)  and 
of  upward  flowing  ion  (UFI)  distributions  has  direct  implications  on  the 
nature  of  auroral  acceleration  processes. 

If  upward  flowing  ions  occur  on  closed  field  lines  they  may  become 
trapped  through  pitch  angle  scattering  and  contribute  to  the  equatorial 
particle  populations.  In  the  absence  of  significant  pitch  angle  scattering 
these  ions  would  appear  in  the  conjugate  hemisphere  as  DFI  distributions. 

Thus  the  occurrence  frequency  and  the  characteristics  of  DFI  distributions 
provide  information  not  only  on  the  interactions  taking  place  during  trans- 
port through  the  equatorial  region  but  also  on  the  injection  of  this  ionos- 
pheric plasma  into  the  equatorial  population. 

Downward  flowing  ion  distributions  in  the  altitude  range  of  '1  R have 
only  recently  been  observed  (Fennell  et  al . , 1978;  Ch lelmet 1 1 o t a 1 . , 1978b). 

At  low  altitudes  (s2000  km)  anisotropical ly  precipitating  ions  have  previously 
been  reported  (Kbr.c  and  Bosqued  , 1971;  Hu  It  q v i s t , 1971,  1978).  Similar  field 
aligned  "source  cone"  structures  of  energetic  ions  (Mauk  and  Mcllw.iin,  1975; 

Borg  c t al . , 1978;  Geiss  et  al . , 1978),  and  bouncing  clusters  of  ions  inferred 
to  be  0 (Me  1 lw.i in , 1976)  have  been  reported  at  geosynchronous  altitude. 

In  this  paper  we  present  the  results  of  a study  of  downward  flowing 
field  aligned  ion  distributions  observed  by  the  S3-3  satellite  in  the  critical 
altitude  range  of  'lR^  where  ions  are  generally  accelerated  upward.  The 
characteristics  of  these  events,  their  relationship  to  magnetospher ic  boundaries 
and  to  the  trapped  particles  arc  investigated  and  the  implications  arc  dis- 
cussed. Two  examples  of  broad  regions  of  downward  flowing  ions  and  enhanced 
plasma  regions  are  described  in  detail,  followed  by  a statistical  study  of  the 
occurrence  of  these  events. 
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EXPERIMENT  DESCRIPTION 


The  measurements  to  be  described  were  made  by  a set  of  3 energetic 
ion  mass  spectrometers  and  A broadband  electron  spectrometers  aboard  the 
S3-3  (1976-63B)  satellite.  The  spacecraft  is  in  an  eliptical  polar  orbit 
with  apogee  at  -8000  km,  perigee  at  250  km,  an  inclination  of  97.5  degrees 
and  a period  of  -3  hours.  Since  the  spin  axis  is  oriented  perpendicular 
to  the  orbital  plane  and  the  sensors  are  mounted  with  their  view  directions 
perpendicular  to  the  spin  axis,  nearly  complete  pitch  angle  scans  are  obtained 
once  every  half  spin  (-9  seconds).  The  three  ion  spectrometers,  each  at  a 
different  energy  setting,  simultaneously  acquire  complete  mass  scans  over 
the  mass/charge  range  from  1 through  -32  AMU/charge  once  every  second. 

Every  16  seconds  the  energy  setting  of  each  spectrometer  is  stepped  through 
one  of  A exponentially  spaced  values,  thus  providing  a 12  point  energy  spectrun 
between  0.5  - 16  keV/charge  in  6A  sec.  The  mass  resolution  is  sufficient 


to  separate  H , He  , He  and  0 ions  at  relative  intensities  of  less  than 


.-2 


10  . The  electron  spectrometers  cover  the  energy  range  from  0.07  < E < 23.5 

keV  in  f our  almost  contiguous  energy  bands  with  a nearly  constant  response 
function.  The  acceptance  angles  of  the  ion  and  electron  spectrometers  are 
approximately  +3°  at  full  width.  More  detailed  descriptions  of  the  instru- 
ment have  been  presented  in  previous  publications  (Johnson  et  al.,  1977; 

Sharp  et  al. . 1978a). 
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OBSERVATIONS 

Approximately  370  orbits  of  data  were  examined  for  downward  flowing 
ion  (DFI)  events.  From  this  data  set  we  have  selected  the  two  most  conspicuous 
DFI  events  for  detailed  discussion.  Both  involve  large  fluxes  of  highly 
anisotropic  ll+  ions  precipitating  over  a wide  latitudinal  region. 

Downward  flowing  ion  event  of  REV  619 

Overview.  This  event  was  observed  near  local  evening  7600  km  above 
the  northern  auroral  zone.  It  occurred  on  September  24,  1976  at  -0020  IT 
during  the  late  recovery  phase  of  a minor  magnetic  storm  that  had  developed 
on  September  18.  This  storm  is  characterized  by  a peak  Kp  of  6 on  September 
18  and  20,  only  moderate  ring  current  injection  with  a peak  list  of  —72-,  , and 
a slow  recovery.  The  sum  of  Kp  was  13+  for  the  day  of  the  observation 
(a  Q-day)  and  20  for  the  precceding  day.  A few  hours  before  this  event 
Dst  showed  a small  negative  excursion  to  -35y  which  is  suggestive  of  a 
substorm  injection. 

A survey  plot  showing  the  raw  detector  counting  rates  from  this  event 
as  a function  of  universal  time  (UT) , altitude  (ALT),  invariant  latitude  (1L\), 
magnetic  local  time  (MLT)  and  L value  is  given  in  Figure  1.  The  top  panel 
indicates  the  energy  step  of  the  3 ion  spectrometers  and  the  period  (16  sec) 
spent  on  each  step.  The  next  two  panels  show  the  logarithms  of  tie  sums  of 
the  counts  from  all  three  spectrometers  from  the  mass  channels  sensitive 
to  ions  with  mass  per  charge  11/Q  - 1 and  16  respectively.  These  values  are 
displayed  once  per  second  corresponding  to  20®  rotation  in  pitch  angle.  The 
latter  is  indicated  in  the  panel  labeled  PITCH,  where  zero  is  the  direction 
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of  precipitation  into  the  ionosphere.  During  this  section  of  the  orbit 
pitch  angles  down  to  -5*  were  sampled.  The  four  bottom  panels  give  the 


| 

l 


logarithm  of  the  count  rates  from  the  four  electron  detectors  sampled  twice 
per  second.  For  a better  characterization  of  the  various  plasma  regimes 


encountered  during  this  section  of  the  satellite  trajectory  we  present  the 
derived  plasma  properties  in  Figure  2.  These  include  the  electron  and 
proton  number  densities,  energy  densities,  average  energies,  and  precipitated 
energy  fluxes  projected  to  ionospheric  altitude.  They  were  calculated  on 
the  assumption  that  the  fluxes  are  concentrated  at  the  logarithmic  centers 
of  the  detector  energy  bands  and  take  into  account  the  measured  pitch  angle 
distribut  ions.  The  periods  of  averaging,  generally  approximately  one  minute, 
and  the  statistical  uncertainties,  unless  smaller  than  the  symbol  size,  are 
indi  catcd  by  the  horizontal  and  vertical  bars  respectively.  We  note  first  in 
Figure  2 a rapid  increase  by  about  an  order  ol  magnitude  in  the  electron 
energy  density  accompanied  by  a decrease  and  a minimum  in  the  average  energy 
at  -71°  ILA.  These  features  are  interpreted  as  the  inner  edge  of  the  plasma 
sheet  (Vasy 1 iunas , 1968,  1970;  F rank , 1971;  SchielJ  and  Frank,  1970).  The 
unusually  low  number  densities  and  high  average  energies  just 
equatorward  of  this  boundary  may  result  from  the  limited  energy  range  of  the 
electron  detectors  (E  > 75  eV) . The  region  of  relatively  high  proton  and 
electron  fluxes  extending  poleward  of  this  boundary  up  to  '80°  ILA  is  identi- 
fied in  Figure  2,  as  the  plasma  sheet.  Further  poleward  the  fluxes  drop  to 
near  the  detector  backgrounds.  The  more  energetic  electrons  (1.6  - 5 koV  and 
7.3  - 23.5  keV)  exhibit  up  and  dovneoming  loss  cones  out  to  -78”  ILA  where 
they  drop  below  the  sensitivity  threshold  of  the  detectors  (see  Figure  l). 
This  boundary  is  indicated  in  Figure  2 as  the  energetic  electron  "trapping 
boundary".  Although  it  is  not  possible  to  unnmiguously  identify  the  boundary 
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of  closed  field  lines  from  particle  signatures  the  presence  of  downcoming 
loss  cones  in  the  pitch  angle  distributions  is  generally  interpreted  as 
Indicative  of  a closed  field  line  geometry.  The  actual  boundary  of  closed 
field  lines  may  be  further  poleward  however. 

For  the  purpose  of  the  following  discussion  we  subdivide  the  data 
period  into  three  principal  regions  based  on  the  proton  plasma  characteristics. 
These  are  designated  as  regions  I,  II  and  III  in  the  upper  margins  of  Figures 
1 and  2.  The  region  designated  I,  extending  inward  of  the  inner  edge  of 
the  plasma  sheet,  exhibits  low  proton  number  densities  and  high  average 
energies.  Region  II  was  chosen  to  cover  the  lower  latitude  section  ('71°  to 
*74®  ILA)  of  the  plasma  sheet,  containing  large  fluxes  of  low  energy  protons. 
Region  III  covers  the  remaining  high  latitude  section. 

A subregion  between  0019:40  and  0021:00  UT  within  Region  II  and  desig- 
nated as  I fa  in  Figures  1 and  2,  is  of  key  interest.  Referring  to  Figure  1, 
we  note  that  the  11  (M/Q  = 1)  fluxes  peak  in  the  downward  direction  (PITCH  = 0”) 
within  this  subregion.  To  quantify  this  peaking  in  the  distributions  we 
define  the  anisotropy  as  the  ratio  of  average  source  cone  (0°  £ u s 30°)  flux 
to  the  average  trapped  (60°  <_  n <_  120°)  flux.  For  the  4 consecutive  satellite 
spins  in  Region  I la  the  anisotropies  ranged  between  2 and  26  with  statistical 
significance  typically  i 3o . The  repetitive  character  of  the  anisotropy  is 
convincing  evidence  that  it  resulted  from  a genuine  pitch  angle  dependence 
and  not  from  temporal  or  spatial  variations  in  the  fluxes.  In  addition  to 
this  extended  region  of  DFI  distributions,  throe  other  isolated  DFI  distri- 
butions occurred  in  the  higher  latitude  portion  of  Region  11  with  anisotropies, 
ranging  between  2.2  and  2.9  and  statistical  significance  above  2o . Between 
the  extended  DFI  event  indicated  by  the  shaded  area  in  Figure  2 and  the 
isolated  DFI  events,  a subregion  of  enhanced  electron  energy  fluxes  and  field 
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aligned  precipitation  of  low  energy  (-1  keV)  electrons  was  encountered.  This 

subregion,  extending  from  approximately  72°  to  72.5*  I LA , is  indicated  by  the 

shaded  area  labeled  UF1  in  Figure  2.  It  also  included  an  upward  flowing 

field  aligned  beam  of  0.5  keV  0 ions  with  peak  flux  of  -9  x 10  (cm  -sec- 

8ter-keV)  1 (Figure  1).  The  significant  decrease  in  the  proton  numb-,  r density 

and  increase  in  the  average  energy  within  this  subregion  resulted  from  a 

reduction  of  the  trapped  low  energy  component  of  the  proton  distribution.  These 

features  are  suggestive  of  downward  electron  and  upward  ion  acceleration. 

Two  additional  multiple  spin  UFI  events  occurred  in  region  III.  The  event 

equatorward  of  the  electron  "trapping  boundary"  was  the  most  intense  with 

8 2 -I 

peak  fluxes  of  -1  x 10  protons  (cm  -sec-ster-keV)  at  0.5  keV.  The  more 

+ + 

poleward  event  involved  both  H and  0 ions  at  energies  of  0.5  - 1.76  keV. 

Angular  and  energy  distributions.  The  average  energy  and  pitch  angle 
distributions  for  regions  I,  II  and  III  are  shown  in  Figures  3 and  4. . The 
differential  energy  spectra  were  obtained  by  averaging  the  fluxes  over  pitch 
angles  from  30°  to  150°  for  the  trapped  component  and  from  0°  to  30°  for  the 
precipitating  component.  The  angular  distributions  were  obtained  by  integra- 
ting the  fluxes  over  the  specific  energy  ranges  indicated  in  Lhe  figures.  We 
see  that  the  proton  fluxes  peaked  at  energies  of  '10  keV  or  above  in  region 
1 (Figure  3a),  and  exhibited  the  typical  upward  and  downward  loss  cones  of 
trapped  fluxes  with  weak  or  no  pitch  angle  diffusion  (Figure  4c).  In  region 
II,  both  the  trapped  and  precipitating  components  of  the  proton  fluxes  increased 
dramatically  at  low  energies  (Figures  3b  and  4d)  resulting  in  an  energy  spec- 
trum peaked  at  *1  keV.  The  angular  distribution  of  these  low  energy  protons 
was  peaked  well  within  the  downward  loss  cone.  This  peak  is  primarily  caused 
by  the  multiple  spin  DFI  event  in  subregion  Ila  since  the  fluxes  were  approxi- 
mately isotropic  in  the  downward  hemisphere  for  the  remaining  section  of 
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region  II  (Figure  5).  The  precipitated  energy  flux  resulting  from  the 
DFI's  in  this  subregion  is  strongly  enhanced  (Figure  2)  and  sufficiently 
Intense  to  excite  visible  proton  auroras.  The  more  energetic,  ring  current- 
like, component  of  the  proton  flux  remained  comparable  to  that  of  region  I 
except  that  the  downward  coming  loss  cone  (0°)  was  filled  in.  At  higher 
latitudes,  in  region  III,  the  fluxes  decreased  over  the  entire  spectrum, 
but  most  dramatically  at  low  energies  (Figure  3c),  and  became  isotropic  in 
the  downward  hemisphere  (Figure  4g,h,i).  The  peak  near  180°  in  Figure  4g 
results  from  the  UFI's  discussed  previously. 

Ion  composition.  The  abundances  of  ion  species  other  than  H+  were 
determined  for  the  low  (0.5  i E < 2.1  keV)  and  high  (5.4  i E < 16  keV)  energy 

parts  of  both  the  trapped  and  precipitating  components.  No  statistically 

2+  + + 

significant  amounts  of  He  , He  or  0 ions  were  detected  in  any  of  the 

three  regions  however.  Note  that  the  upstreaming  0+  event  in  region  II  is 

not  included  in  these  averages  since  only  trapped  and  precipitating  fluxes 

are  considered.  The  limits  to  the  flux  ratios  in  the  low  energy  downward 

-3  2+  + 

streaming  component  of  region  II  were  (3  ± 5)  x 10  for  He  /H  and  (1  i 6) 

x 10”3  for  0+/H+. 

Downward  flowing  ion  event  of  REV  1393 

Overview.  The  second  example  of  downward  flowing  ions  to  be  discussed 

i 

occurred  during  the  main  phase  of  the  December  29,  1976  magnetic  storm  at 
-1051  UT.  The  observations  of  ionospheric  ions  in  the  ring  current  during 
this  time  period  are  discussed  by  Johnson  et  al.  (1977).  This  event  was 
observed  in  the  premidnight  local  time  sector  near  '0300  hours  MLT  at  an 
altitude  of  *5000  km  above  the  southern  auroral  zone.  The  inner  edge  of  the 


J 


bT- 


LMSC/D673078 

plasma  sheet,  although  not  as  clearly  defined  as  in  the  previous  example,  was 
estimated  to  be  at  *56°  ILA  on  the  basis  of  the  decrease  in  the  electron 
energy  density.  The  plasma  sheet  extended  to  *69®  ILA  where  the  more  energe- 
tic electron  and  proton  fluxes  dropped  below  the  sensitivity  threshold  of 
the  instrument.  Both  upward  and  downward  loss  cones  were  observed  in  the 
energetic  electrons  up  to  -67°  ILA,  which  is  indicative  of  closed  field 
lines  in  the  region  equatorvard  of  this  "trapping  boundary".  Downward 
flowing  ions  were  encountered  on  8 consecutive  spins  during  a 2.5  minute 
interval  between  59.5°  to  62.8°  ILA.  They  involved  peak  fluxes  of  up  to 

6.7  x 10  (cm  -sec-ster-keV)  and  pitch  angle  anisotropies  ranging  from 

1.7  to  11.1  with  a statistical  significance  of  typically  i 3o . Although 
significant  fluxes  of  trapped  0+  ions  were  present  in  this  event  only  the 

H+  ions  were  field  aligned  (i.e.  maximum  along  the  magnetic  field  direction). 
As  in  the  previous  case  this  region  of  downward  flowing  ion  distributions 
was  designated  as  region  II.  Regions  I and  III  were  taken  to  extend  from 
57.0®  to  59.5®  and  62.8°  to  65.7°  ILA  respectively.  Thus,  t ho  DFI  event 
was  located  in  the  more  central  section  of  the  plasma  sheet.  A region  of 
downward  auroral  electron  acceleration  approximately  1.5®  wide  occurred  at 
the  electron  "trapping  boundary",  but  the  data  showed  no  evidence  of  UFI 
distributions.  Although  the  plasma  number  and  energy  densities  exceeded 
those  of  the  first  example  (sec  Table  I),  they  were  not  unusual  for  these 
regions  of  the  plasma  sheet  (Vasv I iunas , 1970). 

Energy  and  angular  distributions.  The  differential  energy  spectra  for 
the  three  regions  are  shown  in  Figure  6.  We  note  the  presence  of  a low 
energy  trapped  proton  component  in  regions  I and  II  with  a relatively  flat 
distribution  below  -1  keV.  These  fluxes  decreased  significantly  (by  nearly 
a factor  of  10)  in  region  III.  At  higher  energies  (i 5 keV)  the  ring 
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current  proton  fluxes  gradually  decreased  with  latitude,  but  remained  significant 

in  all  3 regions.  A comparison  of  the  number  densities  and  integral  flux 

ratios  (see  Table  1)  indicates  that  the  low  energy  trapped  protons  were 

enhanced  in  region  II  with  respect  to  the  neighboring  regions.  As  a result 

the  average  energy  minimized  here.  The  angular  distribution  of  these  low 

energy  protons  in  region  II  was  strongly  peaked  in  the  downward  direction 

(180*  pitch  angle  for  this  case  since  it  was  in  the  southern  hemisphere) 

with  a peak  flux  that  was  about  an  order  of  magnitude  higher  than  the  90® 

flux  (Figure  7d).  At  intermediate  and  high  energies  the  angular  distributions 

indicated  onset  of  strong  pitch  angle  scattering  near  the  poleward  edge  of 

region  II.  Thus  the  differential  energy  spectrum  for  the  precipitating 

component  showed  decreased  fluxes  at  higher  energies  and  a broad  peak  at 

about  1 keV  (Figure  6).  Although  the  electron  fluxes  were  generally  isotropic 

outside  the  loss  cones  in  this  region,  strongly  field  aligned  precipitating 

electrons  with  energies  Si  keV  were  observed  on  4 spins.  Examples  of  these 

are  displayed  in  Figure  8.  The  precipitated  energy  flux  from  these  electrons 

2 

projected  to  the  ionosphere  was  *3.5  (erg/cm  -sec)  (on  spin  7)  while  the 

2 

low  energy  protons  precipitated  at  an  average  rate  of  -0.2  (erg/cn  -sec). 

Ion  composition.  Although  H+  was  again  the  major  constituent  of  the 
ion  fluxes,  significant  fluxes  of  0+  and  He+  ions  were  also  observed.  The 
0+  ions  were  more  abundant  at  higher  energies  as  indicated  by  the  differential 
energy  spectrum  for  region  I,  shown  in  Figure  6.  The  0+/H+  flux  ratio  in 
this  higher  energy  component  (5.4  < £ l 16  keV)  decreased  progressively  with 

increasing  latitude  going  from  (1.3  i .1)  x 10  * in  region  I,  to  (7.8  * 1) 

-2  -2 
x 10  in  region  II,  to  (4.1  ♦ .8)  x 10  in  region  III.  At  low  energies 

(0.5  4 E i 2.1  keV)  this  ratio  was  (3.6  ♦ 1.9)  x 10  ^ in  region  I.  A statis- 

+ , + -2 

tically  significant  He  /H  ratio  (2.4  ♦ .6)  x 10  was  obtained  only  for  the 
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high  energy  component  in  region  I.  For  all  other  regions  and  plasma  components 
the  fluxes  of  the  minor  ion  species  were  consistently  near  background,  (e.g. 
within  the  low  energy  downward  flowing  component  of  region  II  these  ratios 
were  (-.9  ± 1.3)  x l(f2  for  He2*/H+,  (.8  + 1)  x 10~2  for  He+/H+  and  (.6  ± .7) 
x 10"2  for  0+/H+). 

Statistical  features  of  downward  flowing  ion  events 

- — — - — — . — — U — - 


To  determine  the  morphology  of  the  DFI  d °nts,  data  in  for  form 
of  survey  plots  of  the  type  shown  in  Figure  1 were  scanned  visually  for  sig- 
natures of  field  aligned  anisotropies.  The  study  was  limited  to  the  high 
latitude  regions  poleward  of  the  radiation  belts  (usually  >60°  I LA)  and 
to  altitudes  above  *2000  km.  It  was  based  on  the  same  data  set  that  was  used 
in  the  statistical  study  of  UFI  distributions  (Chielmetti  et  al. , 1978a), 
covering  the  time  period  from  July  1976  to  February  1977  and  including  -370 
satellite  orbits. 

For  the  purpose  of  this  survey  we  defined  the  pitch  angle  anisotropy 
as  the  ratio  of  flux  in  the  downward  direction  to  the  flux  at  '90°.  In 
general,  a downward  field  aligned  component  must  be  distinguished  from 
statistical  fluctuations  of  the  isotropic  population.  The  estimated  sensi- 
tivity thresholds  for  this  survey  arc  as  follows:  the  pitch  angle  anisotropy 
Bust  exceed  *7.5,  4 and  2.5  if  the  isotropic  fluxes  are  *1,  7 and  20  x 10^ 

i 

2 

(keV/cm  -sec-ster-keV) , respectively.  To  be  accepted  as  a DFI  event  we 
required  that  either  anisotropies  meeting  the  above  criteria  recur  for  at 
least  two  consecutive  spins  (multiple)  or  that  the  anisotropy  exceed  *10  on 
single  spins  (single).  With  these  constraints  only  2 multiple  and  6 single 
spin  DFI  events  were  observed;  corresponding  to  a frequency  of  occurrence 
of  less  than  .05  per  orbit. 
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These  DFI  events  were  observed  at  altitudes  above  -4000  km  on  the  night 
side  (1500  - 0700  MLT).  They  typically  consisted  of  field  aligned  H+  ions 
only,  with  energies  less  than  *3  keV.  Occasionally  more  energetic  protons 
with  energies  up  to  16  keV  and  0+  ions  were  observed.  DFI  events  were 
frequently  located  cquatorward  of  and  at  times  immediately  adjacent  to  UFI 
or  auroral  electron  acceleration  regions.  At  least  one  event  (REV  2210)  was 
discovered  where  DFI  distributions  occurred  on  several  spins  within  a wide 
region  of  UFI  distributions.  In  this  case  the  UFl's  consisted  largely  of 
0+  ions  with  energy  less  than  -4  keV  while  the  DFl's  were  exclusively  H+ 
with  energies  up  to  16  keV. 

The  local  time  region  around  dawn  was  inadequately  sampled  during 
the  time  period  of  the  initial  study.  Therefore,  to  investigate  whether 
or  not  a dawn-dusk  asymmetry  of  the  type  observed  in  the  UFI  events  extisted, 
a separate  study  of  data  acquired  when  the  orbit  of  the  S3-3  satellite  was 
in  the  dawn-dusk  plane  was  made.  For  this  purpose  data  from  May  to  Septem- 
ber 1977  were  added.  A total  of  200  orbits  satisfied  the  condition  of 
passing  within  *2  hours  MLT  of  dusk  and  dawn.  In  this  separate  study  we 
identified  9 single  spin  events  and  1 multiple  spin  event  in  the  dusk  sec- 
tor and  7 single  spin  events  and  no  multiple  spin  event  in  the  dawn  sector. 
Thus,  in  contrast  to  the  UFI  events  (Ghielmetti  et  al.,  1978),  no  statis- 
tically significant  difference  is  indicated  in  the  occurrence  frequency  of 
DFI  events  between  dawn  and  dusk. 

DFI  events  were  often  observed  in  association  with  subregions  of 
enhanced  fluxes  of  trapped  ions  of  similar  energies.  To  determine  the 
frequency  of  occurrence  of  these  subregions,  data  from  the  dusk  portion  of  the 
previously  described  separate  study  were  examined.  For  this  analysis  we 
defined  an  "enhanced  plasma  event"  as  a large  (-10x)  energy  flux  increase  in 


C-14 


LMSC/D673078 


the  soft  protons  (E  s 5 keV)  relative  to  the  harder  (E  i 5 keV)  protons. 

In  addition  we  required  that  the  more  energetic  protons  be  present  over 
latitudinally  extended  regions  both  equatorvard  and  poleward  of  these  sub- 
regions.  Using  this  criterion  any  plasma  enhancements  located  immediately 
adjacent  to  the  low  or  high  latitude  boundary  of  observable  energetic  proton 
fluxes  would  have  been  excluded.  The  study  was  restricted  to  the  dusk  sector 
where  the  inner  boundary  of  the  plasma  sheet  could  be  identified  from  the 

decrease  in  the  electron  energy  fluxes  (see  for  example  REV  619,  Figure  1). 

The  results  of  this  study  are  presented  in  Figure  9.  The  shaded  histogram 
shows  the  probability  of  observing  an  enhanced  plasma  event  within  the  indicated 
latitude  interval  in  the  16-20  MLT  sector  above  6000  km  altitude.  For  com- 
parative purposes  the  latitudinal  distribution  of  UFI  events  averaged  over  the 
two  dusk  local  time  sectors  has  been  included  from  Chielmetti  et  al.  (1978a). 

Enhanced  plasma  events  were  observed  on  approximately  25%  of  the  orbits  and 

most  frequently  between  63°  and  74°  1LA.  Their  low  latitude  boundary  was 
on  the  average  located  - 1 L-unit  poleward  of  the  inner  edge  of  the  plasma 
sheet  and  their  average  width  was  ~1.7L-units.  The  majority  of  DFI  events 
observed  were  found  to  occur  within  these  enhanced  plasma  subregions.  Since 
both  types  of  events  occur  relatively  infrequently  and  both  types  have 
relatively  narrow  latitudinal  extension,  this  correlation  is  significant  and 
further  indicative  of  a causal  relationship  between  the  two  phenomena. 
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DISCUSSION 

In  both  examples  of  downward  flowing  ions  (DFI)  presented  above  the 
distributions  were  field  aligned  in  pitch  angle  and  peaked  in  energy. 

These  features  are  suggestive  of  an  acceleration  preferentially  parallel 
to  the  magnetic  field.  Several  mechanisms  have  been  proposed  to  account 
for  such  distributions  in  the  auroral  electrons  (see  reviews  by  Evans , 

1975;  FAlthammar . 1977).  Fermi  acceleration  has  been  discussed  in  detail 
by  Sharbcr  and  Heikkila  (1972).  This  process  requires  multiple  bounces 
between  hemispheres  and  thus  results  in  both  upward  and  downward  loss  cones 
in  the  ion  and  electron  fluxes.  In  contrast  the  pitch  angle  distributions 
described  above  for  the  protons  (Figures  5 and  7)  and  for  the  field  aligned 
electrons  (Figures  1 and  8)  are  peaked  well  within  the  downward  loss  cone, 
and  thus  are  inconsistent  with  having  resulted  primarily  from  Fermi  accelera- 
tion. 

Alternatively,  parallel  electric  fields  may  give  rise  to  field  aligned 

distributions.  If  the  DFI  are  produced  by  electrostatic  acceleration  of  an 

initially  isotropic  plasma,  the  degree  of  field  alignment  observed  requires 

2 

a relatively  cold  source  population  with  energies  typically  si  x 10  eV 
(near  the  satellite)  and  decreasing  with  increasing  altitude.  This  excludes 
the  hotter  plasma  sheet  population  as  a source.  It  is  not  possible  on  the 
basis  of  the  ion  distributions  to  uniquely  distinguish  between  parallel 
acceleration  above  but  near  the  satellite  and  at  the  conjugate  hemisphere 
however. 

The  detectability  of  a DFI  event  resulting  from  downward  acceleration 
above  the  spacecraft  is  a function  of  both  the  degree  of  field  alignment 
(anisotropy)  and  the  flux  amplitude.  A precise  determination  of  the  effective 
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probing  distance  for  downward  directed  fields  above  the  spacecraft  would 

require  a knowledge  of  many  plasma  parameters  that  are  not  presently  known; 

however,  estimates  of  the  probing  distance  can  be  obtained  on  the  basis  of 

some  simplifying  considerations.  For  a model  of  parallel  electrostatic  and 

scatter-free  acceleration  such  as  discussed  by  Kaufmann  et  al..  (1976), 

the  minimum  plasma  density  required  to  provide  detectable  fluxes  at  the 

satellite  may  be  calculated.  Assuming  a Maxwellian  plasma  with  typical 

-3 

temperatures  of  <10  eV  the  required  densities  are  *1  cm  . In  the  absence 

of  pitch  angle  scattering  and  perpendicular  heating  these  resultant  field 

aligned  fluxes  will  be  contained  in  the  downward  hemisphere  if  the  source 

is  within  -8  R^,  above  the  satellite.  However,  in  view  of  the  observed 

angular  width  of  UFI  distributions  within  the  acceleration  region  (10-20° 

FVHM)  this  simple  model  may  not  be  realistic.  If  one  uses  pitch  angle 

distributions  with  these  greater  widths  immediately  after  acceleration, 

and  takes  into  account  the  sensitivity  levels  discussed  earlier,  the  maximum 

effective  probing  distance  above  the  satellite  reduces  to  a more  typical 

value  of  -3  R^.  The  typical  ambient  plasma  densities  and  temperatures 

measured  on  the  S3-3  satellite  at  altitudes  up  to  8000  km  are  5-50  cm  ^ and 

1-5  eV  respectively  (Moser  et  al . , 1978).  The  ambient  plasma  densities 

measured  during  the  DFI  events  on  REV  619  and  REV  1398  were  -10  cm  ^ and 
2 -3 

-10  cm  respectively  (F.  S.  Mozcr,  private  communication) . Values  of 
-1  cm  ^ are  not  unusual  for  the  regions  at  or  beyond  the  synchronous  orbit 
(Lennartsson  and  Reasoncr,  1978;  Curnett  and  Frank.  1974).  Since  the  am- 
bient plasmas  conditions  arc  generally  adequate  to  produce  detectable  DFI's 
ve  conclude  from  the  low  frequency  of  occurrence  of  these  DFI  events  that 
if  downward  directed  quasistatic  parallel  electric  fields  involving  poten- 
tial drops  * 500  V occur  in  the  altitude  range  from  2000  km  to  -3  Rj.  they 
arc  an  infrequent  phenomenon. 


• C-17 


I 


1 


LMSC/D673078 


The  above  conclusion  implies  a strong  asymmetry  in  the  mech- 
anism for  parallel  ion  acceleration,  favoring  the  upward  direction.  Such 
an  asymmetry  may  indeed  be  expected  if  the  parallel  electric  fields  arc 
primarily  related  to  the  effects  of  magnetic  mirroring  as  discussed  by 
Lennar tsson  (1976).  According  to  other  models  parallel  electric  fields 
may  result  from  current  driven  plasma  instabilities.  Since  the  peak  cur- 
rent densities  of  the  large  scale  upward  and  downward  flowing  currents 
observed  by  the  TRIAD  satellite  (Iijiina  and  Potemra,  1976,  1978)  are  of 
comparable  magnitude,  it  appears  that  the  downward  flowing  currents  may  be 
less  effective  in  driving  such  instabilities  in  this  altitude  range. 

Evidence  that  downward  acceleration  of  positive  particles  may  occur 
at  lower  altitudes  has  been  reported  by  Ri>mr  and  Bosqued  (1971)  and  Hu  1 1 - 
qvist  (1971,  1978).  However,  the  anisotropic  ion  fluxes  observed  by  Hu  1 1- 
qvlst  et  al . (1971,  1978)  exhibit  different  characteristics  from  the  OKI 
events  presented  here  suggesting  a different  process.  Intense  beans  of 
narrowly  collimated  upward  flowing  field  aligned  electrons  have  been 
observed  near  the  magnetic  equator  (lie II wain,  1975)  and  more  recently  at 
the  S3-3  orbit  (Sharp  et  al . , 1978b).  They  have  not  been  observed  simul- 
taneously with  DFI  events  and  appear  to  be  a separate  phenomenon.  At  an 
altitude  of  *1  R^.  Fennel  1 et  al  . (1978)  have  recently  reported  the 
observation  of  a DFI  event  that  was  interpreted  as  consistent  with  a 
local  downward  directed  parallel  electric  field.  However,  the  obser- 
vations of  a region  of  downward  electron  acceleration  within  the  region 
of  DFI  in  REV  619  and  concurrent  with  DFl's  in  REV  1938  do  not  appear  to 
be  consistent  with  a simple  model  of  a quasistatic  downward  directed 
parallel  electric  field  located  above  but  near  the  satellite. 

The  two  examples  of  DFI  events  presented  in  detail  in  this  paper 
occurred  in  the  inner  regions  of  the  plasma  sheet  on  closed  field  lines. 
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Examination  of  the  plasma  parameters  in  Figure  2 and  Table  1 and  of  the 
differential  energy  spectra  in  Figure  3 and  6 reveals  that  both  events  were 
located  within  a wider  region  of  strongly  enhanced  low  energy  (<2  kcV) 
trapped  proton  fluxes  (enhanced  plasma  events).  The  enhancement  was  parti- 
cularly large  in  the  case  of  REV  619  where  the  number  density  increased  by 

2 

factors  of  -8  and  >10  relative  to  the  two  neighboring  regions.  Although 
the  increase  was  relatively  smaller  in  the  example  of  REV  1398,  possibly 
because  of  the  larger  fluxes  of  hot  ring  current  plasma  present  during  the 
magnetic  storm,  it  was  still  significant.  The  more  energetic  (>5  keV)  ring 
current  protons  exhibited  relatively  unstructured  fluxes  throughout  the 
plasma  sheet  and  remained  unaffected  during  these  enhanced  plasma  events. 
Thus  these  subregions  appear  as  high  density  plasma  clouds  of  relatively 
cold  (kT  '.5  keV)  protons  with  limited  latitudinal  extension  (~ 3°  ILA) . 

They  are  embedded  within  the  plasma  sheet  and  within  the  more  extended 
regions  of  hotter  ring  current  fluxes.  The  enhanced  plasma  cloud  ions  and 

concurrent  DFl's  were  observed  at  similar  energies  (Figure  3,  6).  Since 
their  angular  distributions  were  in  both  cases  wider  (see  Figure  5b,  7d) 
than  the  loss  cones,  a fraction  of  the  downstreaming  population  was  trapped. 


In  order  to  investigate  the  relationship  between  the  DFl  and  the 
trapped  components  in  more  detail  vc  subjected  the  anisotropy  and  the  ratio 
of  the  low  to  high  energy  components  of  the  trapped  fluxes  in  region  II  to 
a linear  regression  analysis  on  a spin-by-spin  basis.  In  the  first  example 
of  REV  619  the  correlation  was  not  found  to  be  significant;  however,  this 
data  set  contained  only  4 spins  with  DFl's  and  the  statistical  uncertainties 
were  large.  Using  the  8 spin  periods  in  region  II  and  the  first 
and  last  spin  of  regions  III  and  I respectively  in  the  example  of  REV  1398 
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vc  obtained  a correlation  coefficient  of  0.8.  This  value  is  significant 
•t  the  12  level  suggesting  that  a correlation  did  in  fact  exist  between  the 
field  aligned  precipitating  and  trapped  component  of  the  low  energy  fluxes. 

The  ionic  composition  of  the  fluxes  provides  a further  means  of  dis- 
tinguishing between  the  origins  of  the  different  components.  In  the  example 

of  REV  619  only  protons  were  positively  identified;  the  absence  of  signifi- 
2+ 

cant  fluxes  of  Up  ions  in  the  precipitating  component  is  consistent  with 
an  ionospheric  origin  for  the  plasma.  In  the  event  of  REV  1398  the  high 

energy  component  contained  substantial  amounts  of  0+  and  He+  ions  while  the 
enhanced  low  energy  plasma  was  devoid  of  these  species  and  exhibited  essen- 
tially the  same  composition  within  statistical  limits  as  the  DF1  population. 

The  results  discussed  above  are  thus  suggestive  of  a causal  relation- 
ship between  the  trapped  and  the  field  aligned  precipitating  low  energy 
components  and  are  consistent  with  a common  origin.  Since  the  field  aligned 
component  has  a higher  phase  space  density  it  is  a potential  source  for  the 
enhanced  plasma  cloud.  It  is  therefore  suggested  that  in  these  two  examples 
both  the  DFI's  and  the  enhanced  trapped  plasma  clouds  resulted  from  the 
injection  of  UFl’s  in  the  conjugate  auroral  ionosphere.  The  absence  of 
significant  fluxes  of  0 ions  in  the  DFI  component  may  have  resulted  from 
the  initial  composition  of  the  source  IT  I and/or  from  velocity  dispersion. 
Since  upward  streaming  0+  ions  have  generally  lower  fluxes  and  longer 
bounce  periods  than  the  accompanying  H+  ions,  their  filling  times  become 
considerably  longer  and  transport  processes  correspondingly  more  effective. 

llpwsrd  flowing  ions  arc  known  to  be  a commonly  occurring  phenomenon 
in  the  high  altitude  auroral  ionospheres  (Chiolmcttl  et  nl.,  1978a).  They 
•re  generally  highly  anisotropic  and  involve  peak  fluxes  that  arc  up  to  one 
order  of  magnitude  higher  than  those  of  the  DFI's  in  the  examples  presented 
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(Shelley  et  al«,  1977;  Mlzera  and  Fe  nnell,  1977;  Chielmettl  et  al.  . 1978a). 
Bouncing  clusters  of  Ions  inferred  to  be  0+  (Me II wain.  1976)  and  field 
aligned  energetic  ions  have  been  observed  at  the  geosynchronous  orbit  by 
DeForcst  and  Me  II wain  (1971),  Hauk  and  Me  I Twain  (1975),  and  more  recently 
by  Geir.s  ct  al.  (1978),  and  Bo r g e t a 1 . (1978).  These  observations  suggest 
indeed  that  some  UF1  cross  the  equatorial  plane  prior  to  becoming  isotropic 
and  may  thus  reach  the  conjugate  high  altitude  ionosphere  where  they  would 
appear  as  DFI's  or  enhanced  plasma  events,  depending  on  the  degree  of  pitch 
angle  scattering  during  transit. 

A comparison  between  the  frequency  of  occurrence  of  DFI  and  UFI  events 
indicates  that  the  latter  occurred  approximately  lOx  to  lOOx  more  frequently 
than  the  former  during  the  same  time  period.  Only  a fraction  of  this  differ- 
ence may  be  attributed  to  the  somewhat  lower  sensitivity  for  detection  of 
DFI  dist  ributions.  Thus  we  conclude  that  pitch  angle  scattering  is  usually 
effective  within  1/2  bounce  period  in  isotropizing  the  UFI  distributions, 
and  further  that  the  injection  of  energetic  ions  from  the  auroral  UFI  source 
Into  the  trapped  population  must  be  a common  process. 

The  time  required  for  a UFI  source  to  fill  a flux  tube  isotropically 

to  a specified  flux  level  can  be  estimated  by  comparing  the  total  particle 

contents  of  the  flux  tube  to  the  upgoing  flux.  Assuming  that  the  UFI's  are 

completely  isotropized  within  1/2  bounce  period  and  using  typical  values  for 
8 2 

the  flux  (-1  x lb  /cm  -scc-ster-keV)  at  1 keV  and  typical  widths  of  the 
angular  distributions  (10-20°  FWilM)  we  calculate  filling  times  of  the  order  of  1 
hour  to  reach  the  Isotropic  flux  levels  observed  in  the  regions  of  enhanced 
plasmas  in  the  two  examples  discussed  previously.  However,  considerably 
ahortcr  filling  times  (of  the  order  of  minutes)  would  result  from  the  most 
intense  UFI  events  that  have  been  observed.  These  estimates  do  not  include 
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the  effects  of  losses  from  c.g.  charge  exchange  and  longitudinal  and  ra- 
dial transport;  however,  these  effects  arc  not  expected  to  be  important  in 
cases  of  rapid  flux  tube  filling.  Although  the  actual  durations  of  UFI 
events  are  not  known,  times  of  the  order  of  1 hour  arc  not  unreasonable  con- 
sidering the  duration  of  other  related  phenonena  such  as  auroral  substorms. 

Thus  UFI  fluxes  are  sufficient  to  inject  clouds  of  trapped  ions  on  time  scales 
and  with  densities  comparable  to  the  substorm  associated  equatorial  plasma 
clouds  (DeForest  and  Mcllwain,  1971).  The  proposed  injection  process  could 
result  in  a cloud  of  ions  of  lat itudinally  limited  extension  with  energies 
corresponding  to  those  of  the  UFI's  (typically  a few  keV) . The  process 
would  not  directly  affect  an  equatorial  population  such  as  the  more  energetic 
ring  current  component  that  was  transported  inward  from  regions  deeper  in 
the  plasma  sheet  and  could  thus  result  in  multiple  component  plasmas  of 

the  type  observed.  Enhanced  plasma  events  as  described  above  are  often 
observed  in  the  same  general  regions  as  UFI  events  in  the  dusk  local  time  sec- 
tor, but  peak  at  somewhat  lower  latitudes  as  shown*  in  Figure  9.  The  difference 
in  the  latitudinal  distributions  of  UFI  and  enhanced  plasma  events  may  result 
from  radial  transport  combined  with  the  must  longer  filling  times  at  higher 
L-values.  The  latitudinal  location  of  the  peak  in  the  occurrence  frequency 
of  enhanced  plasma  events  is  also  seen  to  be  in  good  agreement  with  the  quiet 
time  (Kp  1 2)  substorn  injection  boundary  inferred  by  Me  I lwain  (19  74)  for 
the  dusk  region.' 

These  facts  strongly  suggest  that  some  of  the  enhanced  plasma  events 
observed  at  low  altitudes  correspond  to  the  substorm  Injected  plasma  clouds 
in  the  equatorial  plane,  and  that  they  result  at  least  in  part  from  the  in- 
jection of  upward  flowing  ions  from  the  ionosphere.  I'FI  events  are  known  to 
occur  in  a latitudlnally  well  defined  zone  with  maximum  probability  of  occur- 
rence in  the  dusk  hemisphere  (Cl*  le  lmet  t i c t al . , 1978a).  For  a qualita- 
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tive  illustration  we  show  In  Figure  10  a mapping  of  the  local  time-latitude 
distribution  of^UFl  events  into  the  equatorial  plane  using  a simple  dipole 
approximation.  The  substorm  injection  boundary  from  Mcl Iwain  (1974)  has 
been  included  for  comparison.  We  note  that  in  this  projection  the  location 
of  the  most  frequently  occurring  UFI  events  generally  corresponds  to  the  re- 
gion of  the  injection  boundary,  although  the  peak  frequency  of  occurrence  is 
at  somewhat  higher  L values.  Thus  the  auroral  ionosphere  may  be  a signifi- 
cant contributor  to  the  trapped  ions  that  dominate  the  dynamics  of  the 
equatorial  magnetosphere  in  the  vicinity  of  the  geostationary  orbit.  The 
statistical  association  of  DF1  events  with  enhanced  plasma  events  further 
supports  these  suggestions. 

These  results  are  in  agreement  with  the  conclusions  of  Sharp  et  al . . 
(1976a,  1976b)  that  a local  ionospheric  injection  process  acts  over  wide 
range  of  local  times  and  latitudes  during  geomagnetic  storms.  Near  equator- 
ial measurements  by  Frank  (1970),  DoForest  and  Mcl lvain  (1971)  and  Mcl Iwain 
(1974)  have  indicated  plasma  injection  events  well  inside  of  the  geosta- 
tionary orbit  during  magnetically  disturbed  periods.  The  example  of  REV  1398 
corresponds  to  such  an  injection  at  low  L-values  of  -4.  Structured  intense 
fluxes  of  low  energy  ions  have  been  previously  observed  in  the  near  equator- 
ial region  beyond  L * 8 (Frank,  1967)  and  at  low  altitudes  (Frank  and  Acker- 
son,  1972).  These  exhibit  similar  characteristics  to  the  enhanced  plasma 
regions  and  DFI  events  presented  here  and  may  be  related  to  them.  The 
antisunward  flowing  ions  observed  in  the  magnetotail  (Frank  et  al. , 1977; 
Hardy  et  al . , 1977)  may  result  from  UFI's  that  were  injected  in  the  higher 
latitude  regions. 
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SUMMARY  AN!)  CONCLUSIONS 

Enhanced  clouds  of  ions  with  fluxes  and  energies  similar  to  the 
equatorial  plasma  clouds  (DeForest  and  Mcllw.iin,  1971)  are  frequently 
observed  at  altitudes  of  -1  in  the  auroral  zones  by  the  S3-3  sate- 
llite. In  the  dusk  local  time  sectors  these  regions  cup  tquatori ally 
to  the  region  of  the  substorm  injection  boundary.  Downward  flowing 
field  aligned  ions  have  occasionally  been  detected  in  the  auroral  zones. 

They  occurred  preferentially  within  these  enhanced  plasma  regions. 

Two  examples  of  downward  flowing  ion  ( DF 1 ) and  enhanced  plasma  events 
were  described  in  detail.  The  correlations  observed  between  the  two  pop- 
ulations are  suggestive  of  a common  origin  for  both  phenomena  in  the  up- 
ward flowing  ion  (UFI)  events  occurring  in  the  conjugate  hemisphere. 

DFI  events  are  a much  less  frequent  phenomenon  than  UFI  events,  implying 
that  injection  of  upward  flowing  ions  into  the  trapped  population  is  a 
common  and  significant  process.  The  local  tine  and  latitude  distributions 
of  UFI  and  enhanced  plasma  events  and  their  similar  energies  further  sug- 
gest a relationship  to  the  equatorial  plasma  clouds. 

The  results  described  in  this  paper  are  consistent  with  the  direct 
injection  of  spatially  localized  clouds  of  hot  (O.b-5  kcV)  ionospheric 
ions  into  the  plasma  sheet,  from  the  auroral  zone  upward  flowing  ion  source. 
A significant  fraction  of  the  equatorial  substorm  injected  plasma  clouds  may 
thus  result  from  this  ionospheric  source  region.  The  more  energetic  ring 
current  component  ( i 5 kcV)  may  have  originated  in  part  from  injection  of 
UFI  on  higher  latitude  field  lines  and  subsequent  energization  by  inward 
transport. 
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Figure 


Survey  plot  for  the  dusk  section  of  REV  619  on  September  24,  1976 
showing  ion  and  electron  count  rates  as  a function  of  universal  time. 
The  count  scales  are  logarithmic  and  each  tic  mark  on  the  vertical 
scales  represents  one  decade.  The  various  regions  indicated  in  the 
top  margin  are  described  in  the  text. 


Figure  2:  Electron  and  proton  number  density,  energy  density,  average 

energy  and  precipitated  energy  flux  projected  to  the  ionosphere  are 
given  as  a function  of  invariant  latitude  for  REV  619.  Values  are 
averaged  over  periods  indicated  by  horizontal  bars.  Statistical 
errors  are  indicated  by  vertical  bars  when  larger  than  symbol  size. 


Figure  3:  Differential  energy  distributions  for  trapped  (circles)  and  for 

precipitating  (triangles)  H ions  averaged  over  regions  I,  II  and 
111  of  REV  619.  Background  has  been  subtracted. 


Figure  4:  Pitch  angle  distributions  for  low,  intermediate  and  high  energy  H+ 
ions  in  regions  I,  II  and  HI  of  REV  619.  Number  fluxes  were  inte- 
grated over  the  energy  ranges  indicated.  The  dashed  lines  correspond 
to  background  levels.  0®  represents  the  downward  direction. 


Figure  5: 


Figure  6: 


Pitch  angle  distributions  for  low  energy  H ions  in  region  Ila  and 

I 

II  - Ila  of  REV  619.  0°  represents  the  downward  direction. 

Differential  energy  distributions  for  trapped  (circles)  and  prccipi- 

+ ' 

tating  (triangles)  H ions  in  regions  I,  II  and  III  and  for  trapped 
(solid  circles)  0+  ions  in  region  1 of  REV  1398.  Background  has 
been  subtracted. 
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7:  Pitch  angle  distributions  for  low,  intermediate  and  high  energy 
H+  ions  in  regions  I and  1 1 of  REV  1398.  Number  fluxes  were  inte- 
grated over  the  energy  ranges  indicated.  The  dashed  lines  correspond 
to  background  levels.  180#  represents  downward  direction. 

8:  Examples  of  pitch  angle  distributions  of  low  energy  electron  in 
region  II  of  REV  1398. 

9:  Histogram  distribution  of  the  probability  of  observation  of  upward 
flowing  ion  events  and  of  enhanced  plasma  events  in  the  altitude 
range  from  6000-8000  km,  as  a function  of  invariant  latitude. 

The  UFI  distribution  was  averaged  over  the  2 local  time  sectors 
indicated. 

10:  Projection  into  the  equatorial  plane  of  the  regions  of  high  prob- 
ability of  occurrence  of  UFI  events  using  a simple  dipole  approx- 
imation. The  substorm  injection  boundaries  inferred  by  Mcllwain, 
(1974)  are  indicated  by  the  heavy  lines  for  Kp-2  and  5. 
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